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This dissertation describes the space of heterochnic orbits for a class of semihnear 
parabolic equations, focusing primarily on the case where the nonlinearity is a 
second degree polynomial with variable coefficients. Along the way, a new and 
elementary proof of existence and uniqueness of solutions is given. Heterochnic 
orbits are shown to be characterized by a particular functional being finite. A 
novel asymptotic-numeric matching scheme is used to uncover delicate bifurcation 
behavior in the equilibria. The exact nature of this bifurcation behavior leads to 
a demonstration that the equilibria are degenerate critical points in the sense of 
Morse. Finally, the space of heterochnic orbits is shown to have a cell complex 
structure, which is finite dimensional when the number of equilibria is finite. 
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Introduction 
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1.1 The Grand Plan 



This dissertation presents some recent progress towards a rather lofty (and very 
difficult) goal. Specifically, there is great interest in understanding the topology of 
solution spaces of systems of semilinear parabolic equations, 

^ = Di{ui,U2, ...Up) + fi(t,X,Ui,U2, ...Up) 
< ^ = D2{Ui,U2, ...Up) + f2{t,X,Ui,U2, ...Up) (I'l) 
...^ = Dp{Ui, U2, ...Up) + fp{t, X, Ui, U2, ...Up) 

where Ui G C^(]R, C'''"(M")), Di are densely-defined linear elliptic (diffusion) oper- 
ators, and fi satisfy reasonable smoothness conditions. Of course, a major problem 
for anyone interested in (11. ip is existence and uniqueness of short-time solutions! 
Although the existence and uniqueness in general is daunting, many interesting 
and important problems have the form (11.11) . (Fortunately, for some special cases, 
existence and uniqueness can be proven, as is done in Chapter |2l) The applications 
of (II. ip are numerous; for instance: 

• The Navier-Stokes equations, which describe fluid flow, can be put in the 
form of (ll.ip [20J. Understanding the topology of the solution space of the 
Navier-Stokes equations gives insight into the onset of turbulence. This has 
applications to fluid ampliflers (viscous fluid logic gates, with no moving 
parts) [21], in which a particular geometry and set of boundary conditions 
allows several semistable equilibria. The orbits which connect these equilibria 
(the heteroclinic orbits or heteroclines) can involve turbulent flows. As a 
result, understanding the topology of the space of heteroclines for such a 
situation might provide insight into turbulence phenomena. 

• Many chemical reactions are of this form, in particular those describing com- 
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bustion. The precise nature of the ignition of a flame is encoded in the 
topology of the solution space. 

• The combination of the Navier-Stokes equations with combustion equations 
can model turbulent combustion phenomena. Such turbulent, reacting flows 
are important in modeling the inside of internal combustion engines. The 
ignition of a turbulent combustion depends very delicately upon the exact 
nature of the flow, and a topological description for such events is lacking. 

• Related to the Navier-Stokes and chemical reaction equations are nonlinear 
wave equations. Many nonlinear wave equations are of the form (11.11) . and 
traveling waves appear as heteroclines connecting two equilibrium states. 
Traveling waves are often stable, in that perturbations of them tend to "wash 
out" over time. However, there are interesting situations where traveling 
waves are suddenly suppressed as a parameter is changed slightly. This is a 
consequence of an abrupt change in the topology of the solution space. 

• In population biology, ( 11. ip describes a number of competing or cooperating 
species. One can ask about the kinds of bifurcations in stable populations 
when a new species is introduced, or when harvesting patterns are changed. 

• One of the most pressing issues in population biology is that of non-native 
invasive species, which disrupt the ecology of many parts of our planet. In 
agriculture, they cause significant crop losses and threaten our protected 
areas. One of the important problems concerning invasive species is how to 
displace them with minimal ecological impact. Understanding the topology 
of the space of solutions for (11.11) would help find optimal control algorithms 
for eliminating (or limiting) the spread of invasive species. There is a vast 
literature on this subject, going back to Fisher [T3j, Kolmogorov, Petrovski, 
and Piskunov j27j. 



3 



1.2 Specialization to a scalar gradient equation 



Since the general setting of (11.11) is much too difficult to allow any kind of progress 
at this time, we must instead consider more specialized situations. To this end, we 
restrict attention to the case of 

• a scalar equation (p=l in ( 11. ip ). 

• in one spatial dimension, where 

• the diffusion operator is the Laplacian, Di = and 

• the reaction term is a polynomial. 

In other words, consider 

^^^^ = An(t, x) + J2 aiix)u'{t, x) = Au + P{u), (1.2) 

where t G M and x G MJ^, and the a, are bounded and smooth. In this dissertation, 
we consider eternal solutions, those u satisfying (11.21) that lie in C"'^(R, C°'°(]R)) for 
< a < 1. In particular, observe that A : C*^'"(]R) — > C°'°(M) is densely defined 
when a is not an integer. 

This kind of equation provides a simple model for a number of physical phe- 
nomena. First, choosing the right side to be Am — + aiu results in an equation 
which can represent a model of the population of a single species with diffusion and 
a spatially- varying carrying capacity, ai(x). As a second application, this equation 
is a very simple model of combustion. If ai is a positive constant, then the equa- 
tion supports traveling waves. Such traveling waves can model the propagation of 
a fiame through a fuel source. 
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1.3 Discussion of the literature 



Equations of the form (11.21) have been of interest to researchers for quite some 
time. Existence and uniqueness of solutions on short time intervals (on strips 
{—to, to) X R") can been shown using semigroup methods and are entirely standard 
[44] . However, there are obstructions to the existence of eternal solutions. Aside 
from the typical loss of regularity due to solving the backwards heat equation, 
there is also a blow-up phenomenon which can spoil existence in the forward- 
time solution to (II. 2p . Blow-up phenonmena in the forward time Cauchy problem 
(where one does not consider t < 0) have been studied by a number of authors [IB] 
[10] [12] [26] [6] [16] HZ]. More recently, Zhang et a/. ([15] [38] [13]) studied global 
existence for the forward Cauchy problem for 

du 

— = Au + uP - V(x)u 
ot 

for positive u, V. Du and Ma studied a related problem in [9] under more restricted 
conditions on the coefficients but they obtained stronger existence results. In fact, 
they found that all of the solutions which were defined for all t > tended to 
equilibrium solutions. 

Eternal solutions to (II. 2p are rather rare. Most works which describe blow-up 
make the assumption that the solution is positive. Unfortunately, blow-up is much 
more difficult to characterize in the general situation, and understanding exactly 
what kind of initial conditions are responsible for blow-up in the Cauchy problem 
for (11.21) is an important part of the question. 

The boundary value problem that results from taking x G i7 C M" for some 
bounded Q (instead of s G M") has also been discussed extensively in the literature 
j20] [2l] [7]. For the boundary value problem, all bounded forward Cauchy problem 
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solutions tend to limits as \t\ oo, and these limits are equilibrium solutions. 

Almost all of the literature (including this dissertation) describing eternal solu- 
tions to (11.21) is restricted to discussing heteroclines. For unbounded domains and 
certain symmetries among the coefficients a^, one can find traveling waves. Since 
the propagation of waves in nonlinear models is of great interest in applications, 
there is much written on the subject. The general idea is that one makes a change 
of variables {t,x) hh> ^ = x — ct which reduces (11.21) to an ordinary differential equa- 
tion. This ordinary differential equation describes the profile of a traveling wave. 
Powerful topologically-motivated techniques, such as the Leray-Schauder degree, 
can be used to prove existence of wave solutions to (11.21) . Asymptotic methods can 
be used to determine the wave speed c, which is often of interest in applications. 
See |l0] for a very thorough introduction to the subject of traveling waves in (II. 2p . 

1.4 A Morse-theoretic approach 

A somewhat less traditional approach to studying (II. 2p exists. This method at- 
tempts to directly compute topological invariants for the space of heteroclinic orbits 
Ti. of (11.21) . It makes use of the fact that equation (11.21) defines the fiow of the 
gradient of a certain action functional, 

A{f)= I \\\Vfr + j^^r^\x)dx. (1.3) 

It is then evident that along a solution u{t) to (11.21) . A(u{t)) is a monotonic function 
in t. As an immediate consequence, nonconstant t-periodic solutions to (II. 2p do 
not exist. This kind of behavior suggests that a Morse-theoretic framework might 
be helpful. 
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Morse theory is concerned with the computation of homotopy or homology 
groups of a Riemannian or Hilbert manifold M by "exploring" it with a suitable 
scalar function / : M ^ M. The function / is selected to satisfy the Morse-Smale 
(-Floer) conditions, namely 

• a nondegeneracy condition: if x G M is a critical point {df{x) = 0), then the 
Hessian at x is nonsingular, 

• the Morse index, which is the number of negative eigenvalues of the Hessian 
is finite for each critical point, 

• stable and unstable manifolds for the gradient flow of / are transverse (the 
Smale condition), and 

• if M is noncompact, there is a compact isolating neighborhood for each pair 
of critical points under the gradient flow [T7] . 

The function / can be thought of as a special kind of "height function" on M. 
One then examines the topology of sets M* = {a; G M\f{x) < t}, which form a 
cover for M. It is straightforward to show that the homotopy class of M* remains 
constant on t G (^0,^1) when there are no critical points in /^^((to, ^i))- The 
homotopy class of M* changes abruptly, however, when f~^{t) contains a critical 
point. Morse theory describes how this homotopy class changes by the attachment 
of handles to M*. A very readable introduction to Morse theory is [30] . 

There is a dual formulation of Morse theory, which uses Witten's complex 
to compute homology instead of homotopy. This approach is better suited to 
understanding differential equations, as it focuses not on level sets, but rather on 
the flow of 

f = V/(x(t)). (1.4) 
7 



Using this flow, one constructs a chain complex (C=k,9*) in which the Ck are free 
modules generated by the critical points of Morse index k. The boundary maps dk 
are then constructed by the formula 



where n{q,p) is the number of heteroclinic orbits of (11. 4p which connect q to p, 
counted with sign. The surprising thing is that this chain complex computes the 
homology of Ml A thorough, modern treatment of Morse theory can be found in 



1.5 Floer homology 

A similar theory can be made to work even if the Morse index of all critical points 
is infinite. Instead of relying on the critical points to supply an index directly, one 
constructs a "relative" index based on the structure of connecting manifolds. This 
theory was first assembled by Floer for the purpose of understanding the homology 
of the space of orbits for an exact symplectomorphism [15]. More recently, Ghrist 
et al. |T9] have done work on a similar theory for a certain evolution of braids. 
What is crucial to Floer theories is that the manifold of heteroclinic orbits which 
connect a given pair of equilibria (a connecting manifold) is finite dimensional, 
and compact modulo time translation. Suppose that the connecting manifold for 
a given pair of equilibria x,y has dimension fi{x,y). One shows that the following 
two relations hold for this dimension: 



where x, y, z are distinct critical points. This relation allows one to assign indices 
/ to the critical points such that y) = l{x) — I{y). Evidently, / is only defined 




p6Cfc_i 



fj,{x, x) = 0, z) = fi{x, y) + ii{y, z), 



(1.5) 
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uniquely up to an additive constant. However, one can use the index / in place 
of the Morse index to construct a Witten complex in the usual way. However, the 
ambiguity in the definition of / means that the degrees of the resulting complex 
are only defined up to this additive constant. In this way, one obtains a kind of 
ambiguous-degree homology theory, which is called "Floer homology." Alterna- 
tively, one may take the dual approach, and use the finite dimensional connecting 
manifolds to assemble a cell complex structure for the space of heteroclines. The 
attaching maps of this cell complex are evidently related to the boundary maps of 
the Witten complex. 

1.6 How to construct a Floer theory for parabolic equa- 
tions 

For the case of semilinear parabolic equations, the following must be established 
in order to construct a Floer theory: 

1. One must compute the Morse index of each critical point, or more properly, 
show that the Morse index is not well-defined due to degeneracy. 

2. One must show that connecting manifolds are finite dimensional, and that 
they form a cell-complex structure for the space of heteroclinic orbits Ti. 

3. The connecting manifolds must obey the additivity relation (11.51) . 

4. The space of heteroclinic orbits 7i must be compact moduluo time transla- 
tion. 

5. One must construct the boundary maps in the Witten complex, and verify 
that they actually form a chain complex that computes the homology of Ti. 
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This dissertation contains proofs of items 1, 2, and most of 3 for a special case of 
(II. 2p . (In the case of a bounded spatial domain, all but item 5 are standard |32j . 
[24].) Rather than working with equation fll.2p in full generality, the later chapters 
use the following special case: 



where tends to zero as |x| — oo. The resulting questions and techniques have 
obvious generalizations to (11. 2p . though there are many technical obstacles. Higher 
degree polynomial nonlinearities may of course have more than two roots, which 
creates the possibility for more complicated equilibrium structure than we analyze 
here. 

This simpler model still provides insight into applications, as it is still a model 
of the population of a single species, with a spatially-varying carrying capacity, 
(f). Indeed, one easily finds that under certain conditions the behavior of solutions 
to (II. 6p is reminiscent of the growth and (admittedly tenuous) control of invasive 
species It is the control of invasive species that is of most interest, and it is 
also what the structure of the boundary maps reveals. In the example given in 
Chapter [HI there is one more stable equilibrium, and several other less stable ones. 
The more stable equilibrium can be thought of as the situation where an invasive 
species dominates. The task, then, is to try to perturb the system so that it no 
longer is attracted to that equilibrium. An optimal control approach is to perturb 
the system so that it barely crosses the boundary of the stable manifold of the 
the undesired equilibrium, and thereby the invasive species is eventually brought 
under control with minimal disturbance to the rest of the environment. 
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1.7 Outline and prerequisite results 

While Chapters [6] and [7| contain the resuhs of most interest to constructing a 
Floer theory for (11.61) . the other chapters provide a number of results that are 
prerequisites. 

1.7.1 Prerequisites that concern the higher degree case 

Chapters [2] and [3] are devoted to the general equation (11. 2p . The first of these 
provides a new proof of short-time existence and uniqueness for (II. 2p . and as a 
side-effect provides a numerical method for approximating its solutions. While 
existence and uniqueness for (II. 2p is standard [44j, the usual proofs are not suited 
to computation. 

The first novel result for (II. 2p is obtained in Chapter [3l where a decay condition 
on the tti allows us to classify heteroclinic orbits of (ll.2p . In particular, if an eternal 
solution u exists and converges uniformly to equilibria as t — ±oo if and only if u 
has finite energy (supremum of the difference in the action functional (II. 3p over all 
time). Without the decay condition on the Oj, finite energy classifies those eternal 
solutions that connect finite action equilibria. 

The result in Chapter [3] is actually quite important for connecting the Morse- 
theoretic results with the analysis. Morse theory, and in particular Witten's com- 
plex, requires the flow to have a gradient structure. As a result, the space on which 
the flow acts must have an inner product structure, so a natural solution space 
would be C^(M, L^(M)). However, the proofs of the cell-complex structure (Chap- 
ter [7]) do not work with timeslices in L^(M). In particular, one needs this space 
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to have a Banach algebra structure (in which the Laplacian is densely defined), 
so the Holder space C°'°'(M) for < a < 1 is more natural. What follows from 
the results of Chapter [3] is that the space of heteroclines Ti lies in the intersection 
C\R, CO'"(M)) n C\R, L\R)), so in fact there is no difficulty (Corollary [30]). 

1.7.2 Prerequisites that concern the quadratic case 

In the remaining chapters (Chapters H] through [8]), only the special case (11. 6p is 
considered. This is quite sufficient to obtain interesting results about the structure 

of n. 

It is important to understand the collection of equilibrium solutions for (\1.Q\\ . 
which are global solutions to 

= Au-u^ + (j) (1.7) 

on all of M. Like (II. 2p . there are obstructions to global existence in (11.71) |39j . 
Indeed, there are fairly few global solutions to (11.71) . We examine solutions to 
this problem under asymptotic decay conditions for in Chapter HI The solution 
reveals delicate bifurcation behavior in the number of equilbria as is varied. 
Further, the asymptotic behavior is such that all global solutions to (II. 7p have 
finite action (see (11.31) ). 

Since they are rare, it is reassuring to construct an example of heteroclinic 
orbits, which is done in Chapter [51 This example makes specific use of the structure 
of the equilibrium solutions, in particular, their asymptotic decay is crucial. 

In order to construct a Morse theory for (II. 2p . understanding the dimension of 
the stable, center, and unstable manifolds of equilibria is important. In Chapter 
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[6] it is shown that the the center/stable manifold's dimension is typically infinite, 
and later in Chapter [7] it is shown that the unstable manifold has finite dimension. 
Each equilibrium solution is in fact unstable, even if its linearization is stable. This 
implies that each equilibrium is a degenerate critical point. This neatly derails any 
hope of using a standard Morse theory, or even using any of its extensions to infinite 
dimensional dynamical systems [31]. (In Chapter [9l Conjecture 11041 suggests that 
restriction of the fiow to Ti. may correct the degeneracy.) 

The most important result of this work is obtained in Chapter U\ where the 
space of heteroclinic orbits is shown to have a cell-complex structure (with finite 
dimensional cells). The dimension of each cell is determined, under a standing 
assumption of transversality (Conjecture 1951) . From the formula for the dimension 
of the cells, it is clear that an additivity rule like (11.51) will hold. This result 
is further explained by an example in Chapter [HI Finally, in Chapter [9], several 
important future directions are outlined. 
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Chapter 2 

Short-time existence and uniqueness 
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2.1 Introduction 




(This chapter has already been pubhshed as [35].) 



Existence and uniqueness of solutions for fl2.ip under reasonable initial condi- 
tions have been known for some time. For instance, [20j and [44J contain straight- 
forward proofs using semigroup methods. The purpose of this chapter is to show 
how a more elementary proof can be obtained from a sequence of explicitly com- 
puted discrete-time approximations. 

Due to their theoretical and computational stability, implicit iteration schemes 
are often prefered over their easier-to-implement explicit analogues. However, in 
the case of semilinear equations, one can form a hybrid implicit-explicit (IMEX) 
method which offers computational and theoretical benefits. The use of IMEX 
methods for approximating semilinear parabolic equations is well-established [2]. 
Many of the recent works on these methods employ discretizations in both space 
and time. These fully discrete approximations can be computed directly by a 
computer. However, one can obtain a stronger condition for convergence of the 
approximation if only the time dimension is discretized [8]. We show how an even 
stronger condition for convergence is met by the Cauchy problem for 



where G L^(]R") fl L°°(]R"), and how convergence of this method provides an 
elementary proof of existence and uniqueness of solutions. 

The Cauchy problem for (12. ip arises in a variety of settings. Notably, some 
react ion- diffusion equations are of this form [12]. Another application is the special 




(2.1) 



i=0 
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case 



du{x, t) 



Au{x, t) — u^{x, t) + a^i^x), 



dt 



where Oq is a nonzero function of x. This situation corresponds to a spatially- 
dependent logistic equation with a diffusion term, which can be thought of as a 
toy model of population growth with migration. 

Following [8], the approximation to be used is 



j=0 

which is obtained by inverting the linear portion of a discrete version of (12.11) . For 
brevity, we shall call (12.21) the implicit-explicit method. (In the summary paper 
[2], this is called an SBDF method, to distinguish it from other implicit-explicit 
methods.) One can compute the operator (/ — hA)^^ explicitly using Fourier 
transform methods, and obtain a proof of the numerical stability of the iteration 
as a whole. 

2.2 A version of the fundamental inequality 

In order to simplify the algebraic expressions, we make the following definitions. 
Definition 1. Let 



oo 





oo 




(2.3) 



i=0 



and 



oo 




(2.4) 



i=0 



Definition 2. Define the analytic functions 



oo 




(2.5) 



i=0 
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and 



oo 



(2.6) 



i=0 



Since we do not discretize the spatial dimension, we can employ some of the 
theory of ordinary differential equations. We therefore first prove a variant of the 
fundamental (Gronwall) inequality for fl2.ip as is done in [23]. The fundamental 
inequality gives a sufficient condition for approximate solutions to converge. A 
slightly weaker version of Lemma [3] was obtained in Theorem 3.1 of |8j, where the 



existence of solutions was required. 

Lemma 3. Suppose {ui}'fZi is a sequence of piecewise functions Ui : [0,T] 
C2(M") n L^iW) n L°°(M"), such that 

1. there exist A,B>Qso that for each i and t G [0,T], \\ui(t)\\i < A and 
IkiWIloo < B, 

2. for each i and t & [0,^], the series gi{\\ui{t)\\i) and goo{\\ui(t)\\oo) converge, 

3. for each t E [0,T], \\-^Ui{t) — -^(^^(t)) ||oo < and limj^oo Q = 0, and 

4. Mi(0) = Ui{0) for alli>0 

Then for each t G [0,T], {uiit)}"^^ is a Cauchy sequence in L^(R"). 

Proof. Let i,j > be given. Let ri{t) = \\ui(t) — Mj(t)||2 = / {ui(t) — Uj{t))'^dx. 
Notice that the fourth condition in the hypothesis gives ri{0) = 0. 



But, II ■^Ui{t)—F{ui(t)) Woo < is equivalent to the statement that for each t G [0, T] 
and X G M", 




F{ui{x, t)) - ei < u[{x, t) < F{ui{x, t)) + 
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giving 



7]'{t) < 2 J {F{u,{t)) - F{uj{t))) {u,{t) - Uj{t))dx 
+2(e. + e,)l\u.(t)-uMdx 

+2iei + ej)\\uiit)-ujit)\\i 

< 2 j i^iuiit) - Ujit))) iuiit) - Ujit))dx 

+2 J iGimit)) - Giujit))) imit) - Ujit))dx + 2(ei + ej)\\uiit) - Ujit)\\i 

< -2J ||V(i.,(t)-i.,(i))|rdx + 2||G'Mi))-G'KW)||2||i.,(i)-i.,W|b 

+2iei + ej)\\uiit)-u,it)\\i 

< 2||G'(i.,(i)) -G'(x.,(i))||2||i.,W -x.,W||2 + 2(Q + e,)||i.,(i) 



Now also 

\\Giu,it)) - Giu,it))h 



Y,aMit)-u)it)) 



k=0 



< J2\Moo\\u^it)-u';it)l 



k=0 

00 

k=0 



{u^ix, t) — Ujix, tff' dx 



00 „ /k-i \^ 

< J]ll«fcllocA / (M.(a;,t)-M,(a;,t))M dx 

k=0 \ \m=0 / 



fc=0 

00 



fc-1 



Y.n^it)u]-'^-\t) 



m=0 



\Uiit) - Ujit)\\2 



vfc=0 



< g'MWu^it) - uM\2, 
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which allows 



V'it) < 2gUB)\\u,it)-uM\l + ^i^i + ^j)\M)-^M\i- 



V'{t) - 2g'^{B)rj{t) < 2(e, + e,)||n,(t) - u^{t)\\. 



so (recall ri{0) = 0) 



2(e, + e,) / e- 

^0 



■^3'^(^HuJs)-uAs)\Uds 



< 



2(ej + ej) / ||ui(s) - Uj{s)\\ids 



< A{ei + tj)Ate 



Hence as i,j oo, r](t) — > for each t. Thus for each t, is a Cauchy 

sequence in L^(]R"). □ 

Remark 4. Since C2(M")nLi(M'^)nL°°(M") C L2(R") and is complete, Lemma 
|3] gives conditions for existence and uniqueness of a short-time solution to (12.11) . 

Lemma 5. Suppose {ui(t)}'^i is the sequence of functions defined in Lemma\^ 
and that u{t) = limj„,oo Wj(t) in L^(R"). Then 

u'{t,x) = lim u'i(t,x) for almost every x, (2.7) 

i— >oo 

wherever the limit exists. 
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Proof. Notice that since each Ui{t) G L°°(R") and ||Mi(t)||oo < B, the dominated 
convergence theorem allows for each x G M" 

\im u[{t, x)dT = lim / u[{T,x)dT 

= lim {ui(t, x) — Ui{0,x)) 
= u{t, x) — ^(0, x) for almost every x. 

Hence, by differentiating in t, 

u'{t,x) = lim u'i(t,x) for almost every x. 

j— >oo 



2.3 The implicit-explicit approximation 

In this section, we consider the case of a 1-dimensional spatial domain, that is, x G 
M. There is no obstruction to extending any of these results to higher dimensions, 
though it complicates the exposition unnecessarily. 

As is usual, the first task is to define the function spaces to be used. Initial 
conditions will be drawn from a subspace of L^(]R)nL°°(]R), as suggested by Lemma 
[31 and the first four spatial derivatives will be prescribed, for use in Lemma [TOl 

Definition 6. Let 

W = nc\R), 

where we interpret C^(]R) as being the space of bounded functions with four con- 
tinuous bounded derivatives. For the remainder of this chapter, we consider the 
case where each of the coefficients G W. Then let X = {/ G l^|5'i(||/||i) < 
oo and fifoodl/lloo) < oo}. We consider the case where the initial condition is drawn 
from X. 
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An approximate solution given by the implicit-explicit iteration will be the 
piecewise linear interpolation through the iterates computed by (12. 2p . A smoother 
approximation will prove to be unnecessary, as will be shown in Lemma [Til 

Definition 7. Suppose /o and h > are given. Put 

/„+i = (/-/iA)-i(/„ + /iG(/„,)). (2.8) 

The function 

uit) =(^1-(J-- nit)^ ^ nit)^ (2.9) 

where n(t) = , is called the implicit-explicit iteration of size h beginning 
at /o. 

Calculation 8. We explicitly compute the operator (J — hA)~^ using Fourier 
transforms. Suppose 

(/ — hA)u{x) = u{x) — hAu{x) = f{x). 
Taking the Fourier transform (with transformed variable uj) gives 

u^to) + huPu^uj) = f{uj), 
muo] 



The Fourier inversion theorem yields 

/I [ Q^^i^-y) \ 

12^ / TTT^^'j 



Using the method of residues, this can be simplified to give 

u{x) = ((/ - hA)-'f) {x) = ^l fiy)e-\y-^\/^'dy. (2.10) 
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Calculation 9. Bounds on the and L°° operator norms of (/ — hA) ^ are now 
computed. First, let / G L°°(R). Then 



{{I-hA)-'f)i 



x] 



< 

< 
< 



-\y-x\/Vh^y 



2Vh 
Vh Jo 



so ||(/-/iA)-i|U < 1. 



Now, let / G L\R). So then 



ll(/-/^A)-V|| 



< 



—00 
1 



2Vh 



00 

00 fOO 



dx 



\f{y)\e-\y--\/^Hydx 



00 J —00 

00 fOO 



< -= / / e 

V/^ ^-00 Jo 



-\y-x 



'^''^dxdy 



< 



\fiy)\dy=\\fh 



which means ||(/ — hA) ^||i < 1. 



The third condition of Lemma [3] is a control on the slope error of the approxi- 
mation. A bound on this error may be established for the implicit-explicit iteration 
as follows. 

Lemma 10. Suppose fo E X , h > 0. Put f{x,t) = fo{x) + tD{x), where 

(I-hA)-\f, + hGifo))-fo 



D 



Then for every < t < h, 



h 



\\f{t)-F{f{tmoo = o{h). 



{2.U] 
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Proof. Recall every function in X will have bounded partial derivatives up to 
fourth order from Definition [6l 



\\f{t)-F{m)i 



< 



< 



< 



D-{A{fo + tD) + G{fo + tDm^ 

/ oo 

D- ^A{fo + tD) + J2<^^ifo + tDy 



1=0 

oo 



i=0 \j=Q 



i=0 



+ 0(/i) 



(/ - hA) 



-1 



h 



-fo - Afo 



{{I-hA)-'-I)G{fo)\+0{h) 



Now, using the fact that (/ - hA)-^ - I = (I - hA)-^{hA), 

||/'(t)-F(/(t))|U < \\{I-hA)-'Afo-Afo 

+{I-hA)-\hA)G{fo)\\^ + 0{h) 

< ||(/ - hA)-\hA)iAf, + G(/o))||oo + Oih) 

< /i||(/-/iA)-i(AF(/o))||oo + 0(M 

< h\\il - /iA)-ioo||(AF(/o))||oo + Oih) = Oih) 



□ 

Lemma 11. Suppose < hi 0. Let Ui he the implicit- explicit iteration of size 
hi beginning at f^ E X on t E [0, T] . Then provided there exist A, i? > such 
that for each i and t G [0,T], < A and ||uj(t)||oo < B, then the sequence 

converges pointwise to a function in t. The limit function is piecewise 
difjerentiable in t. 
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Proof. Let Ui be the implicit-explicit iteration of size hi. By Lemma [TOl the slope 
error is bounded: 

||<(t)-F(M,(t))|U=0(/i,) = Q. 

Notice that 0. Then, since X C C^(M"), Lemma [3] applies, giving a pointwise 
limit function u{t). Finally, since the slope error uniformly vanishes. Lemma [5] 
implies that the solution is piecewise differentiable. □ 



2.4 "yl priori estimates" for the approximate solutions 



Now we demonstrate that the implicit-explicit method converges for all initial 
conditions in X. Specifically, for each /o G X, there exist A,B>0 such that for 
each i and t G [0,T], ||'Uj(i(:) ||i < A and ||Mi(t)||oo < B, given sufficiently small T. 
We begin by recalling that from Calculation [HI the L°°-norm of (/ — hA)^^ is less 
than one. This means that for the implicit-explicit iteration. 



||/n+l||oo < 
< 

< 

< 
< 



fn + hG{f, 



fn I 



h 



71/ II oo 

OO 



i=0 

oo 



fn 1 1 oo ~)~ ^ ^ ^ 1 1 '^i 1 1 oo 1 1 fn I 



i=0 
oo 



fn 1 1 oo ~l~ ^ ^ ^ 1 1 1 1 oo 1 1 fn \ 
fn 1 1 oo ~l~ hcjoQ ( 1 1 fn 1 1 oo ) 



oo 



Hence the norm of each step of the implicit-explicit iteration will be controlled by 
the behavior of the recursion 



fn+l — fn~\~ hgooifn 



(2.12) 
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for /„, h,a > 0. Since we are only concerned with short-time existence and unique- 
ness, we look specifically at h = T/N and < n < A^, for fixed T > and 

N eN. 

Remark 12. The recursion defined by (12.121) is an Euler solver for 

^=g^{y), with 1/(0) = /o. (2.13) 

at 

This equation is separable, and Qoo is analytic near /o, so there exists a unique solu- 
tion for the initial value problem (12.131) for sufficiently short time. Also, whenever 
yit) > 



the function y(t) is concave up. As a result, the exact solution to (I2.13P provides 
an upper bound for the recursion (12.121) . More precisely, we have the following 
result. 

Lemma 13. Suppose y{0) = /o > m (12.131) . Let T > be given so that y is 
continuous on [0,T], and let N eN. Then for each < n < N , fn < y{T), where 
fn satisfies with h = T/N. 

Proof. Since the right side of (I2.13P is strictly positive, the maximum of y is 
attained at T on any interval [0, T] where y is continuous. Furthermore, since 
y{0) > 0, it follows from Remark [T2] that y is concave up on all of [0,T]. There- 
fore, y is a convex function on [0,T]. Hence Euler's method, (12.121) . will always 
underestimate the true value of y. Another way of stating this is that 

fn < y{nh) < y{T). 

□ 
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Using Lemma [131 the growth of iterates to (]2.12p may be controlled inde- 
pendently of the step size. This provides a uniform bound on the sequence of 
implicit-explicit approximations. 

Lemma 14. Suppose < hi = T/i for i E N. Let Ui he the implicit- explicit 
iteration of size hi beginning at fo E X on t E [0,T]. Then there exists a B > 
such that for each i and t G [0,T], we have \\Ui{t)\\oo < B for sufficiently small 
T > 0. 

Proof. Suppose /j„ is the n-th step of the implicit-explicit iteration of size hi. If 
we let ?/(0) = ll/olloo, Lemma [T3] implies that for any i and any Q <n <i 

\\frn\\oo<y{T) 

for sufficiently small T. Hence by (12.91) and the triangle inequality, ||Mj(t)||oo < B 
for alH and t G [0,T]. □ 

With the bound on the suprema of the approximations, we can obtain a bound 
on the 1-norms. 

Lemma 15. Suppose < hi = T/i for i E N. Let Ui he the implicit- explicit 
iteration of size hi beginning at fo E X on t E [0,T]. Then there exists an A > 
such that for each i and t E [0,T], we have ||'Uj(t)||i < A for sufficiently small 
T > 0. 

Proof. First, notice that Lemma [H] implies that there is a i? > such that for 
each i and t E [0,T], we have ||nj(t)||oo < A for sufficiently small T > 0. Again 
suppose fin is the ra-th step of the implicit-explicit iteration of size hi. Then we 
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compute 



i,n+l II 1 



< 
< 

< 

< 

< 

< 
< 



\fi, 



i + /i.||G(/,„)||i 

oo 

l + hiJ2hkftn\\l 
k=0 
oo „ 

i + h,Y^ / \(^kfL\dx 

k=0 

oo 

1 + /ii^ ||/m||^lafc||oo||/m||l + hi\\ao\\i 
k=l 

oo \ 

l + h,Y,\MooB'-'] +h,\\ao\\, 

+ /ii||ao||i 
{1 + hiC) + hiWaoWi 



k=l 



i{l + ^9ooiB)-^\\ao\\ 



This recurence leads to 



11/ 



in 1 



< 



< 



< 



< 



< 



ra-1 



o||i(l + /iiC)" + /i,||ao||i5^(l + /i.Cr 



m=0 



< ||/o||i(l + /^iC^)" + /i^||ao|| 



(1 + h,cr - 1 

h,C 



olli + ^lkolli) (1 + W-^||ao||i 



0111 



C 



CT 



^llaolli 



o||i + ^||ao||i) + ^ 



C 



Ikolh 



o||i + ^l|ao||ije^^-^||ao||i = A 



Once again, by referring to (12.91) and using the triangle inequality, it follows that 
||Mi(t)||i < S for alH and t G [0,T]. □ 

Theorem 16. Suppose < hi = T/i for i G N. Let ui he the implicit- explicit 
iteration of size hi beginning at f^ E X on t E [0, T] . Then, for sufficiently small 
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T > 0, the sequence converges pointwise to a function in t. The limit 

function is piecewise differentiahle in t. 

Proof. This compiles the results of Lemma [TTl Lemma [T31 and Lemma [TSl □ 

Remark 17. These proofs can be generalized further to handle all equations of 
the form 

where G is as in (12.41) . If the operator L satisfies 

• L : Li(M) nL°°(M) nC°°(M) (M) n (M) is a sectorial linear operator 

• ||(J - hL)-^\\i < 1 and ||(J - hL)-^\\^ < 1, 

then the implicit-explicit iteration 

/„+i = (/-/iL)-i(/„ + /iG(/„)) 
converges whenever / G X. 

Remark 18. Additionally, the techniques can be easily extended to handle the 
initial boundary value problem 

du{x, t) 



dt 



Am(x, t) + ^ ai{x)u\x, t), hi X e K C M", t > 



i=0 

with u{x, t) = v{x, t) a given Lipschitz function along dK x [0, oo), for K compact 
with smooth boundary. In this boundary term appears in the estimate for 

r]'{t) in Lemma [3l which depends on the Lipschitz constant of v. Additionally, 
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in Definition [7l one defines fn+i to be the unique solution to the hnear elhptic 
boundary value problem 

with fn+i{x) = v{x, nh) for x G OK. 
2.5 Conclusions 

The convergence proof for the implicit-explicit method presented here has a num- 
ber of advantages. First of all, like all IMEX methods, each approximation to the 
solution is computed explicitly. As a result, a fully discretized version (as is stan- 
dard in the literature) is easy to program on a computer. Theorem [16] therefore 
assures the convergence of these fully discrete methods. 

However, since the implicit-explicit method presented here is discretized only 
in time, the convergence proof actually shows the existence of a semigroup of so- 
lutions. As a result, the convergence proof forms a bridge between the functional- 
analytic viewpoint of differential equations, namely that of semigroups, and the 
numerical methods used to approximate solutions. While the existence and unique- 
ness of solutions for (12.11) has been known via semigroup methods, the proof pro- 
vided here gives a more elementary explanation of how this occurs. In particular, 
it approximates the semigroup action directly. 
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Chapter 3 

Classification of heteroclines 
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3.1 Introduction 



(This chapter is available on the arXiv as [Mj.) 

In this chapter, the global behavior of smooth solutions to the semilinear 
parabolic equation (11.21) 

^^^^ = Am(x, t) - + ^ ai{xy{t, x) = Au + P{u), (3.1) 

for {t,x) e RxW = M"+i is considered, where > 2 and e L°°(M") are 
smooth with all derivatives of all orders bounded. 

The main result is that solutions to (13.11) which are heteroclinic orbits con- 
necting two sufficiently regular equilibrium solutions of (13.11) are characterized by 
finite energy (Definition [2U1) . That this characterization is necessary at all comes 
from the fact that the spatial domain of (13. ip is unbounded. For bounded spatial 
domains, all bounded global solutions converge to equilibria [21]. The strength of 
our result comes from the fact that the finite energy constraint makes solutions 
behave rather well. Therefore, this result is much sharper than what has typically 
been obtained in the past, and it applies to more complicated nonlinear terms. 

The disadvantage is that in doing so, we cannot treat some of the more compli- 
cated aspects of the dynamics. In particular, travehng wave solutions do not have 
finite energy. Even though a traveling wave will often converge locally to equilibria, 
at least one of those equilibria will not be admissible in our analysis. On the other 
hand, we can exclude traveling waves if we require that all the coefficients Oj decay 
fast enough and only consider one spatial dimension. Then our result establishes 
an equivalence between the heteroclinic orbits and the finite energy solutions. 

For somewhat more restricted nonlinearities, Du and Ma were able to use 
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squeezing methods to obtain similar results to what we obtain here. In partic- 
ular, they also show that certain kinds of solutions approach equilibria [9]. In a 
somewhat different setting, Floer used a finite energy constraint for solutions and 
a regularity constraint on equilibria to characterize heteroclinic orbits of an ellip- 
tic problem [15]. The techniques of Floer were subsequently used by Salamon to 
provide a new characterization of solutions to gradient flows on finite-dimensional 
manifolds [37] . In this chapter, we recast some of Salamon's work into a parabolic 
setting, and of course work within an infinite-dimensional space. 

3.2 Finite energy constraints 

From Chapter [21 we have that solutions to (13. ip exist along strips of the form 
{t, x) G / X R" for sufficiently small t-intervals /. One might hope to extend such 
solutions to all of R"+^, but for certain choices of initial conditions such eternal 
solutions may fail to exist. Fujita's classic paper [18] gives examples of this "blow- 
up" pathology. We will specifically avoid it by considering only eternal solutions to 
(13. ip . By eternal solutions, we mean those which are defined for all R"+^, have one 
continous partial derivative in time, and two continous partial derivatives in space. 
It should be noted that eternal solutions to (13. ip are quite rare: the backwards- 
time Cauchy problem contains a heat operator, and so most solutions will not 
extend to all of R"+^ 

Definition 19. Our analysis of (13. ip will make considerable use of the fact that it 
is a gradient differential equation. Recall that the right side of (13. ip is the L^(R") 
gradient of the action functional (ll.3p . defined for all / G C^(R") by 

^(/) ^ / Jii v/wf - ^ /g rrr/-(.M.. 
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It is then evident that along a solution u{t) to ( 13. ip . 
dA{u{t)) 



dt 



dA\u(t) 



dt) 



VA{u{t)) 
Au + P{u), 



du 
'dt 
du 
'dt 



du 
'dt 



> 0, 



so A(u{t)) is a monotone function. As an immediate consequence, nonconstant 
t-periodic solutions to (13. ip do not exist. 



Definition 20. The energy functional is the following quantity defined on the 
space of functions R""*"^ — > R with one continuous partial derivative in the first 
variable [t), and two continuous partial derviatives in the rest (x): 



E(u) 



1 



du 



dt 



+ \Au + P(u)\ dxdt. 



(3.2) 



Calculation 21. Suppose u is in the domain of definition for the energy functional, 
then 

2 



Eiu) 



— oo 
oo 



du 



— oo 
oo 



dt 

du 
'dt 

du 
'dt 

du 
'dt 
du 
di 



lA'u + P{u)\ dx dt 



Au - Piu) \ + 2— (Am + Piu)) dx dt 
dt 



Au-P{u)] dxdt+ [ (^,VA{u{t)))dt 



— oo 
oo 



dt 



d 

dt 



A{u{t))dt 



(^-Au- P{u) ] dxdt + A{u{T)) 



T=-oo 
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This calculation shows that finite energy solutions to ( 13. ip minimize the energy 
functional. If a solution to (13.11) is a heteroclinic connection between two equilibria, 
then the energy functional measures the difference between the values of the action 
functional evaluated at the two equilibria. The main result of this chapter is the 
converse, so that finite energy characterizes the solutions which connect equilibria. 

Remark 22. Finite energy solutions to (13.11) are even more rare than eternal 
solutions. However, the set of finite energy solutions is not entirely vacuous, as 
will be shown in Chapter [51 

It is well-known that when equations like (13.11) exhibit the correct symmetry, 
they can support traveling wave solutions [12]. A typical traveling wave solution 
u has a symmetry like u{t,x) = U{x — ct) for some c G M. As a result, it is 
immediate that traveling waves will have infinite energy. On the other hand, they 
also evidently connect equilibria. As a result. Calculation [2T] shows that a necessary 
condition for traveling waves is that there exists at least one equilibrium whose 
action is infinite. In this chapter, we will consider only equilibria with finite action, 
and solutions with finite energy. As a result, we will not be working with traveling 
waves. 

3.3 Convergence to equilibria 

In this section, we show that finite energy solutions tend to equilibria as \t\ — > oo. 
In doing this, we follow Floer in [15] which leads us through an essentially standard 
parabolic bootstrapping argument. 

Lemma 23. Let U (^W andue W^'P{U) satisfy \\D^u\\^ < C < oo for < j < k 
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(in particular, u is bounded). If P{u) = Xlili '^i'^' with ai E L°°{U) then there 
exists a C such that \\P{u)\\k^p < C'\\u\\k,p- 

Proof. First, using the definition of tlie Sobolev norm, 

k k N 



j=0 j=0 i=l 



CLiZL lip. 



Now \D^aiU^\ < Pij{u, Du, D^u) is a polynomial in j variables with constant 
coefficients, which has no constant term. (It has constant coefficients because the 
derivatives of the are bounded.) Additionally, 

i/p 



1/p 

< ||Z^'"m||1, ( / |D%r ) < C'WD^ul 



pi 



so by collecting terms, 

fc TV k 

\\P{u)\\k,p < 5^5^ ll/^WlIp < 5^A,||D%||p < C'\\u\\k,p. 



j=0 i=l j=0 

□ 



The following result is a parabolic bootstrapping argument that does most of 
the work. In it, we follow Floer in [15], replacing "elliptic" with "parabolic" as 
necessary. 

Lemma 24. If u is a finite energy solution to (13. ip with ||-D''m||l°°{(-oo,oo)xV) ^ 
C < oo for < j < k with k > 1 on each compact V C M", then each of 
limt^±oo u{t,x) exists, and converges with k of its first derivatives uniformly on 
compact subsets o/R". Further, the limits are equilibrium solutions to (13. ip . 
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Proof. Define Umit,x) = u(t + m,x) for m = 0,1,2.... Suppose U C 



pn+l 



IS 



a bounded open set and K C U is compact. Let /? be a bump function whose 
support is in U and takes the value 1 on K. We take p > 1 such that kp > n + 1. 
Then we can consider Um G iy^'^(?7) (recall that u and its first k derivatives of u 
are bounded on the closure of U), and we have 

ll^^mll VK'=+i.P(X) < ||/5Mm||H^fc+i.P(C/)- 

Then using the standard parabolic regularity for the heat operator, 

d 



dt 



Let P'{u) = —u^ + Yld=i diU^i noting carefully that we have left out the term. 
The usual product rule, and a little work, as suggested in [37] yields the following 
sequence of inequalities 



2||Wm||H/'='P(C/) 



/3 ( - - A ) + (3P'{um) - l3P'{ur 



vy'=.p(i7) 



< Ci 

+C2||Mm||Ty'='P((7) 

+ C^\I3P' 

< Ci||/9ao||vF'='P(i/) + C'3||'?im||vy'='P(i/)! 



where the last inequality is a consequence of Lemma [221 By the hypotheses on 
u and oq, this implies that there is a finite bound on llMmll vk'^+^'PCx)) which is 
independent of m. Now by our choice of p, the general Sobolev inequalities imply 
that ||Mm||c'=+i-("+i)/p(A') is imiformly bounded. By choosing p large enough, there 
is a subsequence {vm'} C {um} such that Vm' and its first k derivatives converge 
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uniformly on say to v. For any T > 0, we observe 



dv 



dt 



dx dt 



lim 

m'— »oo 



lim 



T 
m'+T 





dVm' 


/ 


dt 





du 


/ 


dt 



dx dt 



dx dt = 0, 



where the last equality is by the finite energy condition. Hence |^| = almost 
everywhere, which implies that v is an equilibrium and that limj^oo u{t, x) = v{x). 
Similar reasoning works for t —oo. □ 



Now we would like to relax the bounds on u and its derivatives, by showing 
that they are in fact consequences of the finite energy condition. 

Lemma 25. Suppose that either n = 1 (one spatial dimension) or N is odd, then 
we have the following. If u is a finite energy solution to (13. ip . then the the limits 
limt^-too x) exist uniformly on compact subsets, and additionally, 

• u is hounded, 

• the derivatives Du are hounded, 

• and therefore the limits are continuous equilihrium solutions. 



Proof. Note that since 



du 



we have that for any e > 0, 



1 



lim 



T+e 



dt 



du 



+ I Am + P(u)\ dxdt < oo. 



dt 



+ \ Au + P{u)fdxdt = 0, 



whence lim^^oo | ^ | = for almost all x. So this gives that the limit is an equilib- 
rium almost everywhere. Of course, this argument works for t ^ —oo. 
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Now in the case of being odd, a comparison principle shows that solutions 
to (13. ip are always bounded. So we need to consider the case with N even. In 
that case, a comparison principle on (13. ip shows that u is bounded from above. 
On the other hand, if is even we have assumed that n = 1 in this case, and 
it follows from an easy ODE phase-plane argument that unbounded equilibria are 
bounded from below. (Here we have used that the coefficients are bounded.) As 
a result, we must conclude that if a solution to (13. ip tends to any equilibrium, that 
equilibrium (and hence u also) must be bounded. 

Now observe that |^| — > as t — > oo on almost all of any compact K cM."', 
and that | ^ | < a < oo for some finite a on {{t, x)\t = 0, x G K} by the smoothness 
of u. By the compactness of K, this means that if 11% 11 r^^u ^ = oo, there 
must be a (t*,x*) such that \im(^t,x)^{t*,x*) = oo. This contradicts smoothness 
of M, so we conclude |^| is bounded on the strip (—00,00) x K. On the other 

hand, the finite energy condition also implies that for each w e M", 

2 



lim 

s—*oo 



whence we must conclude that lim. 



/•oo /• 


du 


/ — 00 J K+sv 


dt 



dx dt = 0, 



du{t,x+sv) 



at 



for almost every t G M and 



X E K. Thus the smoothness of u implies that is bounded on all of R"+^. 

Next, note that since ||^| and u are both bounded, then so is Am. (Use the 
boundedness of the coefficients of P.) Taken together, this implies that all the 
spatial first derivatives of u are also bounded. 



As a result, we have on a bounded equicontinuous family of functions, so As- 
coli's theorem implies that they (after extracting a suitable subsequence) converge 
uniformly on compact subsets of K to a continuous limit. □ 

Corollary 26. Suppose that either n = 1 or N is odd. An eternal solution u to 



38 



(13. ip has finite energy if and only if each of the following hold: 



• each ofU±{x) = lim^^^-too u{t, x) exists and converges with its first derivatives 
uniformly on compact subsets o/M", 

• U± are bounded, continuous equilibrium solutions to (13.11) . 

• and either \ A{U+) - A{U^)\ < oo orU+ = U^. 

Remark 27. If we consider the more limited case of (11. 6p . then the asymptotic 
decay rate for equihbria indicates that all equilibria have finite action (Corollary 
HH]). In this case, Corollary [261 characterizes a// heteroclinic orbits, not just those 
whose action difference is finite. 

Theorem 28. Suppose that n = 1 and that all equilibria have finite action. If u 
is a finite energy solution then it converges uniformly to equilibria as \t\ —>■ oo. 

Proof. Suppose that u tends to equilibrium solutions f± as t ^ ±oo. Suppose 
that this convergence is not uniform, so that there exists an e > such for each T, 
there is a |t| > |T| with either ||M(t) — /_ ||oo > e or ||u(t) — /+||oo > e- We therefore 
postulate the existence of a pair of sequences {tj}, {xj} such that \tj\ oo and 
mm{\u(tj,Xj) — /_(xj)|, \u{tj,Xj) — f^{xj)\} > e for all j. We assume that for each 
j, Xj is chosen so that min{\u{tj,Xj) — f^{xj)\, \u{tj,Xj) — f+{xj)\} is maximized. 
Notice that since u —>■ f± uniformly on compact subsets, we must have \xj\ oo. 

To simplify the discussion, we find an i? > such that for all \y\ > e and 

|a;| > R, 



We assume that \xj\ > R for all j. This condition ensures that the leading nonlinear 
coefficient of (13.11) dominates. 




(3.3) 



i=0 
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We discern three cases, which we can consider without loss of generahty after 
extracting a suitable subsequence of {{tj,Xj)}. In each of the cases, we shall 
perform a coordinate transformation so that the equilibrium to which u converges 
is the zero function. In particular, we start the sums at 1 rather than 0. 

1. Suppose tj +00 and u{tj,Xj) > e > 0. Since Xj is chosen at a maximum 
of u{tj) for each j, we have that ^ Hh^^j) < q i^y ^i^e maximum principle. As 
a result, 

d ^"^ 

-g^u{tj, Xj) = 'o^'^i^ji^j) ~ {^ji^j) + ^ ] ai(xj)u (tj,Xj) 

i=l 

N-l 

< —u^{tj, Xj) + ai{xj)u^(tj, Xj) 

i=l 

< -y, by (m- 



Therefore we conclude that ||M(t)||oo —>■ 



2. Suppose tj ^ —oo and u(tj,Xj) > e > 0. Consider the time-reversed version 
of fl3.ip . namely 



1=1 

The comparison principle works in reverse for this equation! Suppose that 
v(t, x) = U (t) is a spatially constant solution to ( 13. 4p with U(tj) = v(tj,Xj) = 
u{tj,Xj) > for some j. Then, shortly thereafter, ||M(t)||oo > U{t), since 

du{tj,Xj) 



dt 



> f/^(t,)-$^a.(x,)f/^(t,) 



i=l 

> V(t,)>o. 



On the other hand, this rate of growth indicates that u blows up in finite 
time. This contradicts the fact that u is an eternal solution. 



40 



3. Suppose u{tj,Xj) < — e < and that tj 



-oo or to 



-oo. If is 



odd, then this case can be covered by the previous ones, mutatis mutandis. 
Therefore, we assume is even. We assume that the hmit as t ^ +00 of u 

_ _ . . q2. 

< 00 



is the zero function. From Lemma [251 we have a constant A ■ 
which is independent of t. Thus for each tj, we have an upper bound for 
u{tj) which looks hke 



f A 

Uj{x) = min < /ij, u{tj, Xj) + — (x — Xj)' 



(3.5) 



for some fij > 0. In particular, note that fij ^ by the previous cases. 

We show that the forward Cauchy problem for fl3.ll) started with Uj as an 
initial condition blows up for sufficiently large j. By the comparison principle, 
this implies that u cannot be an eternal solution, which is a contradiction. 
This can be shown using the method of Fujita, which we briefly sketch here. 

Apply the coordinate transformation w = u — ^ for some fi > fij > 0. 
Therefore, the initial condition can be made entirely negative, and by the 
previous cases, the solution stays negative for arbitrarily long future time (by 
taking j large). (Notice that it may not remain negative for all future time 
in the case where tj —>■ —00.) This transformation changes (13. ip into 



ow d^w M i sr^ 



w 



i=l 



k=0 



(-1 


(A 







/i 



N~k 



i=k+l 



J 



As is usual for the Fujita method, we choose a solution v to ^ — 
particular, fix T > and choose e > to define 

1 



In 



v/47r(T -s + e) 



4(T-s + e) ^j)^ 
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Then we define J{s) = J Ve{s,x)w{s — tj,x)dx and compute 



dJ 
ds 



< 



dve dw 
—w + v^—dx 

OS OS 



/ v^w^dx + / ttiV^w^dx 
1=1 ^ 



N-l 
A;=0 




N-k 



i=k+l 



i—k 



dx 



N-l „ 

i=i 



ajVfdx 



N-l 
k=0 



\w\ 



V 



N-l 



i=k+l 



I (^i II ooA'' 



i—k 



J 



where we have used Lemma [2^ to bound w, and we have used the assumption 
that is even in the last step. Since Oj decays to zero, the second term can 
be made arbitrarily small for an arbitrarily large s by taking j large as well 
(for fixed T and e). (The second term may eventually grow larger.) The last 
term is a constant, independent of T, e and can be made arbitrarily small by 
taking j large. By (13.51) . for sufficiently large j, J(0) < 0, and for larger j, 
J(0) becomes more negative. Therefore, for a certain T and sufficiently large 
j, J(s) tends to — oo for some < s < T. However, this contradicts the fact 
that w is bounded. 



□ 

Corollary 29. Suppose that n = 1 and that all equilibria have finite action. If u 
is a finite energy solution then it converges uniformly to equilibria as \x\ oo. 



Proof. Really, the only thing that must be noticed is that Theorem [28] shows that 
there is uniform convergence in the time direction. For a given e, there is a T > 
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such that \u{t,x) — f±{x)\ < e for all \t\ > T. However, this means that for 
t G [— T, T], this does not hold. However, [— T, T] is compact, and the proof of 
Lemma [21] indicates that there is uniform convergence to equilibria as ||x|| oo 
on compact subsets. □ 

Corollary 30. The above Corollary implies that the asymptotic spatial behavior 
of heteroclinic orbits is determined entirely by the asymptotic spatial behavior of 
equilibria. In particular, in Chapter^ it is shown that the equilibria for the case 

of (USD, 

dt dx"^ 

with cj) decaying to zero, all lie in L-^(M). As a result, each timeslice of a heterocline 
lies in L^iM.) as well. 



3.4 Discussion 

The point of employing the bootstrapping argument of Lemma [21] is only to ex- 
tract uniform convergence of the derivatives of the solution. As can be seen from 
the proof of Lemma [211 such regularity arguments are unneeded to obtain good 
convergence of the solution only. 

While Corollary [2H] is probably true for all spatial dimensions, the proof given 
here cannot be generalized to higher dimensions. In particular, Veron in j39] 
shows that in the case of P{u) = —u^, there are solutions to the equilibrium 
equation Am — = which are unbounded below and bounded above when the 
spatial dimension is greater than one. This breaks the proof of Lemma [25] that 
the limiting equilibria of finite energy solutions are bounded for N even, since the 
proof requires exactly the opposite. 
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On the other hand, the case of P{u) = —u\u\'^^^ + J2iLo^ is considerably 
easier than what we have considered here. In particular, all solutions to (13.11) are 
then bounded. In that case, the proof of Lemma 1251 works for all spatial dimensions. 
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Chapter 4 

Equilibrium analysis 
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4.1 Introduction 




(This chapter is available on the arXiv as [33], and has been accepted for publica- 

tion in Ergodic Theory and Dynamical Systems.) 

Since the dynamics of solutions to the semilinear parabolic equation (11 ■2p de- 
pend strongly on the equilibrium solutions, it is important to understand the num- 
ber and structure of equilibrium solutions. As will be shown, this is a somewhat 
ill-defined and rather delicate goal. Therefore, to fix ideas and techniques, we shall 
focus on the specific case of the equilibria of (11.61) 



encounter have obvious generalizations to the more general equation. Therefore, 
we are faced with the task of analyzing a nonlinear ordinary differential equation, 
and finding its global solutions. Additionally, the asymptotic properties of such 
solutions will be crucial in Chapters [3l [5], and [71 

Finding global solutions to nonlinear ordinary differential equations on an infi- 
nite interval can be rather difficult. Numerical approximations can be particularly 
misleading, because they examine only a finite-dimensional portion of the infinite- 
dimensional space in which solutions lie. Additionally, the conditions for global 
existence can be rather delicate, which a numerical solver may have difficulty rig- 
orously checking. In situations where there is well-defined asymptotic behavior 
for global solutions, it is possible to exploit the asymptotic information to answer 
questions about global existence and uniqueness of solutions directly. Additionally, 
more detailed information may be provided by using the asymptotic behavior to 
install artificial boundary conditions for use in a numerical solver. The numerical 





where tends to zero as \x 



oo. The resulting questions and techniques we 
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solver can then be used on the remaining (bounded) interval with boundary con- 
ditions that match the numerical approximation to an asymptotic expansion valid 
on the rest of the solution interval. 



In this chapter, we consider the behavior of global solutions satisfying the 
equilibrium equation for (14.11) . namely 



In particular, we wish to know how many solutions there are for a given (p. (There 
may be uncountably many solutions, as in the case where = const > 0.) This 
problem depends rather strongly on the asymptotic behavior of solutions to (14.21) 
as |s| — > oo, so it is useful to study instead the pair of initial value problems 



where G C°°(R). The sets supply the initial conditions for which solutions 
exist to ( 14. 3p for all x > and to f l4.4p for all x < 0, respectively. Solutions to ( 14. 2p 
will occur exactly when Z {~\Z' is nonempty. Indeed, the theorem on existence and 
uniqueness for ODE gives a bijection between points m Z {~\ Z' and solutions to 
(14. 2 p [29] . Since (14. 3 p and (14.41) are related by reflection across x = 0, it is sufficient 
to study (14.30 only. 

Due to the asympotic behavior of solutions to (14. 3p . the methods we employ 
here will be most effective in the specific cases where is nonnegative and mono- 
tonically decreasing to zero. (We denote the space of smooth functions that decay 



= f"{x) - f{x) + 0(x), for all x G M. 



(4.2) 




(4.3) 



and 




(4.4) 
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to zero as C^(]R).) The decay condition on (p allows the differential operators 
in (14.21) through (14. 4p to be examined with a perturbative approach as x becomes 
large, and makes sense if one is looking for smooth solutions in L^(]R) with bounded 
derivatives. 

When is strictly negative, it happens that no solutions exist to (14. 3 p for all 
X > 0. The monotonicity restriction on (p provides some technical simplifications 
and sharpens the results that we obtain. This leads us to restrict to a class of 
functions that captures this monotonicity restriction but allows some flexibility, 
which we shall call the M-shaped functions. 

It is unlikely that we will be able to solve (14. 3 p explicitly for arbitrary 0, so 
one might think that numerical approximations might be helpful. However, most 
numerical approximations will not be able to count the number of global solutions 
accurately. For instance, finite-difference methods are typically only useful for 
finding solutions valid on finite intervals of M. However, one cannot easily infer 
a solution's behavior for large values of |x| when it is only known on a finite 
interval. In particular, global solutions to (14. 3p must tend to zero (Theorem [38l) . 
All other solutions fail to exist for all of M. Worse, the space of initial conditions 
which give rise to global solutions is at best a 1-dimensional submanifold of the 
2-dimensional space of initial conditions (Theorem [551) . Therefore, a typical finite- 
difference solution that appears to tend to zero may in fact not, and as a result 
fails to be a solution over all a; > 0. 

Because of this failure, we need to understand the asymptotic behavior of so- 
lutions to (14. 3 p as we take x oo. Equivalently, since G C^(M), this means 
that we should examine solutions with small. The driving motivation for this 
discussion is that solutions to = f"{x) — p{x) + (j){x) for small behave much like 
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solutions to = f"{x) — /^(x). In the latter case, we can completely characterize 
the solutions which exist on intervals like [xq, oo). 

In Section 14.21 we review what is known about the much simpler case where 
is a constant. Of course, then (14. 3 p is autonomous, and the results are standard. 
In Section 14.31 we establish the existence of solutions which are asymptotic to 
zero. Some of these solutions are computed explicitly using perturbation methods 
in Section 14.41 where low order approximations are used to gather qualitative 
information about the initial condition sets Z and Z'. In Sections 14. 51 and [4. 61 these 
qualitative observations are made precise. Section 14.71 applies these observations 
about Z and Z' to give existence and uniqueness results for (14. 2p . Finally, in 
Section 14.81 we use the information gathered about Z and Z' to provide artificial 
boundary conditions to a numerical solver on a bounded interval, which sharpens 
the results from Section 221 We exhibit the numerical results for a typical family 
of (j), showing bifurcations in the global solutions to (14. 2p . 

4.2 Review of behavior of solutions to = f"{x) — P{x) + P 
It will be helpful to review the behavior of 



where P is a constant, since varying can be viewed as a perturbation on the case 
= P. In particular, we need to compute some estimates for later use. We 
shall typically take P > 0, as there do not exist solutions for all a; if P < 0. 

Lemma 31. Suppose f is a solution to the initial value problem (14. 5 p with f{0) > 




(4.5) 
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Phase plot of 0=f'(x)-f^(x)+9 




-10 -8 -6 -4 -2 2 4 6 8 10 



f(x) 

Figure 4.1: The phase plot of /" — + 9 = 0. Bounded solutions live in a 
small region, the rest are unbounded. 

sfP and /'(O) > 0. Then there does not exist an upper bound on f{x), when x > 0. 
Additionally, if P <0, there does not exist an upper bound on f{x). 

Proof. Observe that for / > ^/P or if P < 

r = f-p> 0. 

Hence, since /'(O) > 0, and /' is monotonic increasing, f{x) is monotonic increasing 
at an increasing rate. Thus it must be unbounded from above. □ 

Definition 32. The differential equation (14.51) comes from a Hamiltonian, namely 

H{fJ') = lf-lf-fP+lp'/'. 

Definition 33. A useful tool in the study of smooth dynamical systems is the 
funnel. Suppose $ is a local flow on a manifold M. A funnel F is a set such that 
if X G -F, then G F for all t > 0. A funnel F with an oriented, piecewise 
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boundary is characterized by having the vector field -^i^x being inward-pointing 
for all X G dF. 

Lemma 34. Suppose f is a solution to the equation (14.51) on M. All bounded 
solutions lie in the funnel 

M = {{f, f)\H{f, f) >Oandf< ^/P}. (4.6) 

Any solution which includes a point outside the closure of M is unbounded, either 
for X > or X < 0. (Note that M is the teardrop- shaped region in Figure ) 



Proof. • M is a bounded set. Notice that H{f, 0) > H{f, /'), or in other words 
within M, 

< 1 ;3 _ 1^,2 _fp^ 2^3/2 < 1^3 _fp^ 2^3/2 

Elementary calculus reveals that this inequality establishes a lower bound on 
/, namely that 

- V3P <f<VP (4.7) 

On the other hand, 



1/1 < ^JIp'/' + Ip - VP < \j\p'" (4.8) 
immediately establishes a bound on /'. 

M is a funnel, from which solutions neither enter nor leave. This is immediate 
from the fact that H is the Hamiltonian, and the definition of M simply says 
that i/(/, /') > 0. This suffices since solutions to (14. 5 p are tangent to level 
curves of H. 



If (/(0),/'(0)) i M then / is unbounded. Evidently if /(O) > VP and 
/'(O) > 0, then Lemma [3l] applies to give that / is unbounded. For the 
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remainder, discern two cases. First, suppose /(O) > \fP and /'(O) < 0. 
Evidently, i/(/(0), /'(O)) = i/(/(0), -/'(O)), so it's just a matter of verifying 
that a solution curve transports our solution to the first quadrant. But this 
is immediately clear from the formula for 

/' = ±^^/3-2/P-2i/(/(0),/'(0)), 

which gives /' = ±/'(0) when / = /(O). The other case is when 
H{f{0),f'{0)) < 0. Then we show that there is a point {^/P,g) on the 
same solution curve, and then Lemma [311 applies. So we try to satisfy 

- 1/ - P=^/2 + = i/(/(0),/'(0))<0 

g' = -2i7(/(0),/'(0))>0, 

which clearly has a solution in g. 

□ 

Lemma 35. If f is a solution to (14. 5 p with /(O) > -\/P, and /'(O) > then there 
exists a C such that Imix^c f{x) = oo. 

Proof. From Lemma [HH we have that / is unbounded, and goes to +oo. Using 
the Hamiltonian, we can solve for 

£ = ±y^/3 - 2fP - 2H{f{0), /'(O)), 

or viewing / as the independent variable, 

dx 1 

±y^|/3-2/P-2if(/(0),/'(0)) 




52 



as / becomes large. Solving this asymptotic differential equation is easy, and leads 
to 

3 

which has an asymptote at a; = C. □ 



4.3 Existence of asymptotic solutions for G C^(M) 

The first collection of results we obtain will make the assumption that tends to 
zero. From this, a number of useful asymptotic results follow. Working in the phase 
plane will be useful for understanding fl4.3p . Of course (14. 3 p is not autonomous, 
but by adding an additional variable, it becomes so. 

Definition 36. We think of fl4.3p as a vector field on M^, defined by the formula 



(4.9) 



P - 

V ' J 

Notice that the first coordinate of an integral curve for this vector field solves (14. 3p . 

Definition 37. Define H{fJ',x) = - \f'^ - f(t){x) + l(j)^/^{x). Notice that 
for constant = P, this reduces to a Hamiltonian for (14.50 . 

Theorem 38. Suppose f is a solution to the problem (14. 3 p where (j) G C(j^(]R). /// 
does not tend to zero as x oo, then there exists a z such that lim^_»2 f{x) = oo. 
Stated another way, if f solves (14.30 for all x > 0, then limj._,oo f{x) = 0. 
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Proof. If / does not tend to zero, this means that there is an i? > such that for 
each xq > 0, there is an x > Xq so that |/(x)| > R. But since tends to zero as 
X — s> oo, for any P > we can find an xi > such that for all x > Xi, \(f){x)\ < P. 
Choose such a P so that the set M in Lemma [M] associated to (14. 5p is contained 
entirely within the strip —R < f < R. We can do this since the set M is bounded, 
and its radius decreases with decreasing P, as shown in (14 .7^ and (14.81) . But this 
means that there is an X2 > Xi such that |/(a;2)| > R- 

Construct the following regions (See Figure : 

/ = {(/,/',x)|/>Pand /'<0}, 

J/ = {(/,/',x)|/>Pand /'>0}, 
///={(/,/',x)|/<-P}, 

and 

IV = > and / > -P and - ^f' - fP+'^P'^' < 

if /< Vp)}. 

The following statements hold: 

• Region I is an antifunnel. Along f = R and /' = 0, solutions must exit. 
Once a solution exits Region /, it cannot reenter. Also, because / > y/P, 
/" = — > — P > 0, solutions must exit Region I in finite x. 

• Region II is a funnel. Along f = R and /' = 0, solutions enter. Now 
/" = — > — P > and /' > 0, so solutions will increase at an 
increasing rate and so, they are unbounded. 

• Solutions remain in Region /// for only finite x, after which they must enter 
Region IV. This occurs since / < —\fP < 0, and so /' always increases. 
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Figure 4.2: The Regions /, //, ///, and IV of Theorem 



Note that for /' < 0, solutions will enter Region III along / = —R, and for 
/' > 0, solutions exit along / = —R. 

Region IV is a funnel. Solutions enter along / = —R and along /' = (note 
that I/I > y/P in both cases). Along the curve boundary of Region IV, we 
have that 



so that solutions enter. 



\ 







V 1 I 



/'(P-0)<O, 



Now suppose {f{x2),f'{x2),X2) G /. After finite x, say at x = x'2, the solution 
through that point must exit Region /, never to return. Then, there is an X3 > x'2 
such that 1/(0:3)1 > R. So this solution has either (/(xs), /'(xs), X3) G // or G ///. 
The former gives the conclusion we want, so consider the latter case. The solution 
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will only remain in Region /// for finite x, after which it enters Region IV, say at 
X = X3. Then there is an X4 > Xg such that |/(a;4)| > R. Now the only possible 
location for (/(X4), /'( ) to be is within Region IL since it must also remain 
in Region IV. As a result, the solution is unbounded by an easy extension of 
Lemma [31] As x becomes large, (p tends to zero, so the solution will be asymptotic 
to an unbounded solution of = /" — But Lemma [3l] above assures us that 
such a solution is unbounded from above, and Lemma [35] gives that it has an 
asymptote. Hence, our solution must blow up at a finite x. □ 

This result indicates that solutions to fl4.3p which exist for all x > are rather 
rare. Those which exist for all x > must tend to zero, and it seems difficult to 
"pin them down." We now apply topological methods, similar to those employed in 
[23] , to "capture" the solutions we seek. The methods we use are due to Wazewski 

m- 

We begin by extending the usual definition of a fiow slightly to the case of a 
manifold with boundary. 

Definition 39. Suppose M is a manifold with boundary. A flow domain J is a 
subset of R X M such that if x G M then = pr^(J fl M x {x}) is an interval 
containing 0, and if x is in the interior of M then is in the interior of Jx- 
{pTi : M X M — i> R is projection onto the first factor) 

Definition 40. A (smooth) flow is a smooth map $ from a fiow domain J to a 
manifold with boundary M, satisfying 

• $(0, x) = X for all X G M and 

• ^{ti + t2,x) = $(^2, x)) whenever both sides are well-defined. 
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Additionally, we assume that flows are maximal in the sense that they cannot be 
written as a restriction of a map from a larger flow domain which satisfles the 
above axioms. We call the curve '■ Jx ^ M deflned by ^x{t) — ^{t,x) the 
integral curve through x for $. 

Definition 41. Suppose $ : J ^ M is a flow on M and x e dM. Then the flow 
at X is said to be inward-going (or simply inward) if is an interval of the form 
[0, a) or [0, a] for some < a < oo. Likewise, the flow at x is outward-going if Jx 
is of the form (a, 0] or [a, 0] for — oo < a < 0. 

Theorem 42. (Wazewski's antifunnel theorem) Suppose ^ : J ^ M is a flow on 
M and that {A, B} forms a partition of the boundary of M such that the flow of 
$ is inward along A and outward along B. If every integral curve of $ intersects 
B in finite time (ie. Jx is hounded for each x), then A is diffeomorphic to B. 

Proof. For each x e ^, Jx = [0,ta;], where tx is the time which the integral curve 
through X intersects B. (We have that ^{tx, x) is outward-going, since Jx is closed, 
so it is in B.) 

Using this, we can deflne a map F : A ^ B hj F{x) = ^{tx,x). F takes A 
smoothly and injectively into B. The smoothness follows from the smoothness of $ 
and that dM is a smooth submanifold. To see the injectivity, suppose F{x) — F{y) 
for some x,y e A, so ^{tx,x) = ^{ty,y). Without loss of generality, suppose 
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< tx < ty. Then we have that 

Fix) = F{y) 
^-tx,F{x)) = ^-tx,F{y)) 
^{-tx,^{tx,x)) = (!>{-tx,^ty,y)) 
^U-U,x) = ^ty-tx,y) 
X = ^{ty-tx,y). 

But the flow is inward at x, so it is also inward at ^{ty — tx, y). This means that 
{ty — tx — e,y) ^ J for every e > 0. But this contradicts the fact that {ty,y) G J 
unless we have ty <tx- As a result, ty — tx, so x — y. 

In just the same way as for F, we construct a map G : B ^ A so that G takes 
B smoothly and injectively into A. Namely, we suppose Jy — [sy, 0] for some Sy, 
and put G{y) = ${sy,y). Notice that by maximality, if there were to be an a; G A 
such that F{x) — y, Sy — —tx- 

Now we claim that G is the inverse of F. We have that 
{GoF){x) = ^{sF(x),F{x)) 

= ^{sF{x),^tx,x)) 

= ^{sF{x)+tx,x) 

= ^{-tx+tx,x) ^X, 



where we employ the remark about Sy above. □ 

Remark 43. We can extend the Antifunnel theorem to a topological space X 
on which a flow $ : J ^ X acts in the obvious way. In that case, there is no 
reasonable deflnition of the boundary of X. However, the notion of inward- and 
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Figure 4.3: Schematic of the region Ri, showing the boundary partition A 
and B. 

outward-going points still makes sense. If we let A be the set of inward-going 
points and B be the set of outward-going points in X, then the conclusion is that 
A is homeomorphic to B. 

Now we employ the Antifunnel theorem to deduce the existence of a bounded 
solution to = f" — p + (p for x > xq for some xq > 0. 

Theorem 44. Suppose < (f){x) < K for all x > Xq for some Xq > and 
< K < oo, and that there exists an Xi > Xq such that for all x > Xi, > 0. 
Then the region Ri given by Ri = {{f,f',x)\H{f,f',x) > 0,x > 0, / < ^y(f){x)} 
contains a bounded solution to = f"{x) — /^(x) + (p{x), which exists for all x 
greater than some nonnegative X2- 

Proof. Without loss of generality, we may take xq = 0, because otherwise solu- 
tions must exit the portions of Ri in {(/, f',x)\x < xq} since the x-component of 



59 



V{f, /', x) is equal to 1. 



If 0(0) > 0, partition the boundary of -Ri into two pieces: A = {{f, /', x)\x = 0} 
and B = {{fJ',x)\H{fJ',x) = 0} (See Figure iJ]). The flow of V is evidently 
inward along A. As for B, notice that VH is an inward-pointing vector field normal 
to B. We compute 



VH-V 



p - m 



/ + v^)0'(x), 



J 



which has the same sign as 0'(x) when / < ^/^(a;) in R^. Finally, we must deal 



with the case where / = \f^{x) G B. But in this case, /' = from the equation for 



so we see that V{^f^{x), 0, x) = (0, 0, 1)-^, so the flow is inward when (f)'{x) < 
and outward when (f)'{x) > 0. This means that the portion of the boundary of Ri 
on which the flow is outward is a disjoint union of annuli. On the other hand, the 
portion of the boundary of -Ri on which the flow is inward is the disjoint union of 
a disk (namely i?i fl {x = 0}) and some annuli. 

We now consider the case of 0(0) = 0, in which case the set A above is just a 
point. Assume without loss of generality that (j){x) is strictly positive for all x > 0, 
so we let xi = 0. Let 



x 



inf {x G (0, oo)|0'(x) = 0} or 



oo if (j)'(x) > for all x > 0. 



In this case, the set {(/,/', x)|0 < x < x'} fl dRi, is a contractible (it may be a 
point if oscillates rapidly as x ^ 0), connected component of the inflow portion 
of the boundary of -Ri. It is obvious that the remainder of the inflow portion of 
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the boundary is homeomorphic to a disjoint union of annuli, since is smooth and 
strictly positive. 

We can apply the Antifunnel theorem to conclude that there is a solution which 
does not intersect either the inflow or outflow portions of the boundary. There is 
a lower bound on the x-coordinate of such a solution, since the x-component of 
V{f,f',x) is equal to 1, and the Region Ri lies within the half-space a; > 0. 
Therefore, there must exist a solution which enters Ri, and remains inside the 
interior of Ri for all larger x. That such a solution is bounded follows from the fact 
that each constant x cross section of Ri has a radius bounded by the inequalities 
fliTfj) and (gSl), and the fact that < K < oo. □ 



4.4 Asymptotic series solution 

Theorem HH ensures the existence of solutions to = /" — + for x sufficiently 
large. However, it does not give any description of the initial condition set Z which 
leads to such solutions, nor does it give a description of the maximal intervals 
of existence. Fortunately, it is relatively easy to construct an asymptotic series 
for solutions to (14.31) . which will provide a partial answer to this concern. In 
doing so, we essentially follow standard procedure, as outlined in [22], for example. 
However, our case is better than the standard situation, because under relatively 
mild restrictions this series converges to a true solution. 

We begin by supposing that our solution has the form 

oo 

/ = 5Z/., (4.10) 

k=0 

where we temporarily assume fk+i -C fk and /o ^ 0, as x — »• +oo. (This assump- 
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tion will be verified in Lemma l45l ) Substituting (I4.10p into (14. 3p . we get 



fc=0 



fk fmfl 



k~in 



+ 



m=0 



= io'-/o +(/r-2/o/i + 0) + $^ 



k=2 



k-1 



f'k ~ 2/o/fc — ^ fmfi 



k—m 



m=l 



We solve this equation by setting different orders to zero. Namely, 



fe-i 



= f^-fl 

= fk~ '^fofk — fmfk-m- 
m=l 

The equation for /o is integrable, and therefore easy to solve. (There are two 
families of solutions for /q. We select the nontrivial one, because the other one 
simply results in f{x) ~ — (j){s)ds dt.) The equations for are linear and 

can be solved by a reduction of order. Thus formally, the solutions are 



f — 6 

^0 ~ (x-d)^ 
h- ' 



{x-d)3 



K+nt-drj; 



6 r°° 0(^) 



ds dt 



(4.11) 



= - /^(^ - df ir ^^^^^S^ds dt, 
for d, K constants. Notice that these constants parametrize the set of initial con- 
ditions Z. 



Lemma 45. Suppose f{x) = ^^=o fk{x) where the fk are given by (14.111) . If there 
exists an M > 0, an R> 0, and an a > 5 such that 
M 



X < 



for all \x - d\ > R > 0, 



{x — d) 

then f{x) is bounded above by the power series 

1 °° 4 ^ 



k=0 



X — d 



(4.12) 



(4.13) 
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Proof. Wc proceed by induction, and begin by showing that the /i term is appro- 
priately bounded: 



< 



< 



{x - df 
1 



{x - df 

1 



/X nc 



{s - df 
M 



{x — dy 



K + 



t (s-d) 

M 



ds dt 
ds dt 



3+a 



(2 + a) (5 -a){x- d^-^ 



Now since \x — d\> R and a > 5, we have that 



< 



\x — d\^ 
\x-d\^' 



\K\ + 



M 



(2 + q;)|5 - a\R'^- 



with 



A = \K\ + 



M 



{2 + a){a - b)R'^- 



(4.14) 



For the induction hypothesis, we assume that < i^J^j^+j with ^Ij > and 
for alH < /c — 1. We have that 



k-i k-i , , 



^ ^ fmfk—m ^ ^ ] 



m=l 



\x - \x - d\ 



2+k-m 



m=l 



SO by the same calculation as for /i, we obtain 



A A, 



(A; + 6)(A;-1) 
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Hence we should take 



A A, 
{k + 6){k-l) 



(4.15) 



Hence we have that 



oo ^ oo 

|/(a;)|<5^|A(a;)|<^-— 5^/1 



k=0 



k=0 



1 



X — d 



Lemma 46. The power series given by 



oo , 

It — 



fe=o 



X — d\'^^ 



with Aq,Ai > given, and 



Efe-l A A A A 

m=l ^m-'^k—m Z_/m=l ^m-^k- 



{k + 6){k-l) A;2 + 5A;-6 
converges for \x — d\ > R if Ai < 8R. 



□ 



Proof. We show that under the conditions given, the series passes the usual ratio 
test. That is, we wish to show 

Ak+i 



lim 



< R. 



Proceed by induction. Take as the base case, k — 1: by the formula for A^, 

A\ A2 A\ Ai „ 
8 ' 8^1 8 - 
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Then for the induction step, 



^k+l _ X]m=l 

^ ~ Akik^ + 7k) 

Ak{k^ + 7k) 
Mk^+7k) 

- Ak{k^ + 7k) 
RAkik"^ + 5k - 6) + AiAk 

- + 7k) 
R{k'^ + 5k-6) + Ai 

- + 7k) 
+ 5k-2) 

- (P + 7k) 

< R, 



since Ai < 8R. Thus 



< R for all k, so the power series converges. □ 



Lemma HHl provides conditions for the convergence of the bounding series found 
in Lemma HSl Hence we have actually proven the following: 

Theorem 47. Suppose f{x) = YlkLo f^i^) where the fk are given by (14. lip . // 
there exists an M > 0, an R> 0, ana > 5 such that (I4.12p holds, and furthermore 

M <8{a + 2){a-5)R°-\ (4.16) 

then the series for f{x) converges for all x such that \x — d\ > R. 

Proof. Combining Lemmas HS] and HHl we find that the key condition is that Ai < 
8R, which by substitution into (14.140 yields 

' ' {a + 2){a- 5)i?°-5 " 
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M(d) for (|)=(x^-0. 12)6"* a=10 




-10 -8 -6 -4 -2 2 4 6 8 10 

d 



Figure 4.4: A typical M{d) function 

But in order to have \K\ > 0, this gives 

0<8ii ^' 



(a + 2)(a -5)i?"-5' 

which leads immediately to the condition stated. □ 

Corollary 48. As an immediate consequence of Theorem\Jl\ the action functional 
A given by (11.31) evaluated at each equilibrium is finite provided the conditions 
(KWf and KWf on hold. 

Remark 49. It is worth noting that if the spatial dimension n > 1 ( fl4.2p is now an 
elliptic partial differential equation), then the asymptotic decay rate will typically 
be slower than that of the series solution given here. As a result, Corollary H51 will 
not hold for higher spatial dimensions. Indeed, whether anything like Corollary HHl 
holds in higher spatial dimensions is an open question. 

Example 50. It is important to notice that the M defined above in Lemma HS] 
can depend crucially upon the value of d and the shape of the curve 0(x). For the 
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1=3erief c-^iiverges^ □=3erie3 inFi;' diverge c=D.1£, a=1D 




Figure 4.5: Series convergence test, for = (x^ — 0.12)e ^^Z^: white = 
series converges, black = series may diverge 

case of 4>{x) = {x^ — c)e~^^/^, a typical plot of M{d) is shown in Figure l474l It 
should be noted that for various values of c, the M{d) function is numerically very 
similar. 

This also means that the condition (I4.16P defines a somewhat complicated 
region over which parameters d, K and R yield convergent series solutions. An 
example with our given function is shown in Figure 1431 Thus it appears that 
our series solution converges if one goes out far enough, and specifies small enough 
initial conditions. 

Remark 51. The convergence of the series solution is controlled by the conver- 
gence of a well-behaved power series. It follows that as the function becomes 
smaller, fewer terms in the series are needed to accurately approximate the solu- 
tion. Indeed, each term in the series solution is asymptotically smaller than the 
previous one. Thus, we can gain some qualitative information from the leading 
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two terms of the series, which are 
6 1 



fix) 



{x — dY {x — dY 
Taking a derivative by x gives 

^ ^ {x- df (x - dY 



/X fC 
{t - df 

/X po 
{t - df 



[s - dy 



-ds dt 



{s - dy 



-ds dt 



+ 



{x - df 



[s - dy 



-ds. 



On the other hand, using the standard expansion for (a + 6)'^/^, one obtains 



X) ~ 



6^ 



+ 



3-6^ 



{x - dY (x - d) 



K+ / {t-dY 



{s - dy 



:ds dt 



63/2 (a;_rf)3 3.62 



+ 



K+ / {t-dY 



{x-dY 2-63/2 {x-dY 

K+ / {t-dY 



{s - dy 



1/2 



-ds dt 



63/2 3 . 1^ 

+ 



{x-dY &'^{x~dY 



t {s-dy 



:ds dt 



which leads to 

_y|;3/^ + (._,)3£_*i£)_,. (4,7) 

Notice that this equation depends only on d, not K. So from this we should expect 
that the initial data for solutions to be confined to a thin region in the plane x = 0. 
This will be confirmed in Theorem [55] 

Additionally, the relation /' = — a/2/3/3/2 holds exactly for the bounded solu- 
tions of = /" - /2. Indeed, in that case, the set Z is {(/, f')\3f'^ = 2/3, /' < 0}. 
So (14.171) indicates that the presence of 7^ will deflect the set Z largely in the 
/' direction. This is exactly what we show in Section 14.61 
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4.5 Restriction to nonnegative and monotonically de- 
creasing 

We now examine what stronger results can be obtained by requiring > 
and (f)'{x) < for all x > 0. This can be expected to provide stronger results, in 
particular because the region Ri employed in Theorem HH acquires a simpler inflow 
and outflow structure on the boundary, and in particular, solutions will exist for 
all X > 0. A collection of four results indicate that all bounded solutions to (14.31) 
lie within a narrow region. 

Lemma 52. Suppose (f){x) > and < for all x > 0. Then the region given 
by Ri = {(/, /', x)\H{f, /', x) > 0, X > 0, / < a/ (f){x)} contains a bounded solution 

to (531). 

Proof. Following the proof of Theorem HH we partition the boundary of -Ri into 
two pieces: A = {(/,/', x)|x = 0} and B = {{f,f',x)\H{f,f',x) = 0}, noting that 
the flow of V is inward along A. Reviewing the computation in Theorem HH the 
flow is outward along all of B. 

Now we employ the Antifunnel theorem, noting that while A is simply- 
connected, B is not. Hence they cannot be homeomorphic, and so there must 
be a solution that remains inside Ri (which evidently starts on A). But the flrst 
coordinate of such an integral curve must obviously be bounded, since the x cross- 
sections of Ri form a decreasing sequence of sets, ordered by inclusion, and the 
cross-section for x = is a bounded set. □ 

Lemma 53. Suppose 0(x) > and (f)'{x) < for all x > 0. Then the region 
given by R^ = {(/, f, x)\H{f, f, x) < 0, |/3 - > 0, x > 0, /' < 0} contains a 
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bounded solution to ( 14. 3p . 



Proof. Partition the boundary of R2 into two pieces: 



A = {if, 0)1/ < 0} U {(/, x)\Hif, /', a;) = 0, / < /' < 0}, 



and 

B = {{f, /', x)\H{f, x) = 0, / > /' < 0}U 

-^r = 0,/'<0}. 

By the calculation in Theorem the flow along A is inward-going. Additionally, 
the flow along the first connected component of B is outward-going. Finally, we 



put S{f, /', x) = — |/'^ and observe that VS* is an inward pointing normal 



If - if" 
vector field to B. We compute 



vs-v 



-f 

\' J 

rm < 



v 



so the flow along this component of B is outward-going. As a result, we can apply 
the Antifunnel theorem, noting that A is connected, while B is not. Therefore, 
there exists a solution to (14.31) that remains in R2. Note that there is a lower bound 
on the x-coordinate of this solution, since the x-component of V{f, f, x) is equal 
to 1, and the Region R2 lies within the half-space x > 0. So this solution must 
enter R2 through A, and then never intersect B. Additionally, notice that such a 
solution will have /' < and / > 0, so it must be bounded. □ 

Lemma 54. Suppose 0(x) > and (f)'{x) < for all x > 0. The complement of 
the set A = RiU R2 consists of solutions which are unbounded, and blow up in 
finite x. 
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Figure 4.6: The Regions /, //, and /// of Lemma [Ml 

Proof. Let the complement of the set A be called (7, namely C = {(/,/', > 
0} — A. Now the calculations in Lemmas f52l and f53l show that (7 is a funnel, in that 
the flow through the entire boundary of C is inward. If (p does not tend to zero, 
then the argument in the proof of Theorem [38] completes the proof, as there is a 
tubular neighborhood about {/ = /' = 0} with strictly positive radius in which 
solutions in C cannot remain. So without loss of generality, we assume ^ ^ 0. 

Deflne the Region / by 

/ = {(/, x)\f> and (/' > or H{f, /', x) > 0) } . 

There are two bounding faces of Region /, along which the flow is inward. The 
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first is 5*1 = / — a/ 4>{x) = 0, along which 

/ 

VS,-VifJ\x) = 



1 




f 



f 



\ 



> 0. 



The second was computed already in the proof of Theorem HH Notice that /" = 
P — (t>{x) > in Region /, so /(x) is concave- up, so solutions which enter Region 
/ are unbounded. Using similar reasoning to that of Theorem [38l such solutions 
blow up in finite x. 

Now suppose we have a point (a,a',Xo) G C with a' < 0. We claim that for 
some Xi > xq, the integral curve through this point will cross the /' = plane. To 
see this, construct Region // by 

// = <! (/, /', - - ^a' + < and /' < o| n C. 



Note that 



V(^f -^r)-n/,/',a;) = /V(x)<0, 



so the flow is inward along Region II except along /' = (along which it is 
outward). Also note that Region // excludes a tubular neighborhood of the line 
/ = /' = with strictly positive radius. As a result of this, the integral curve 
through ( a, a' , Xo ) proceeds at least as far as to allow / < — a/0(x), at which point, 
a finite amount of distance in x takes it to /' = 0. 

So at that point, the integral curve has entered Region ///, say at x = Xi, 
where 

/// = {(/, /', x)\H{f, f, xi) < and / < and /' > 0}. 
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The flow is evidently inward along /' = and the curved portion by previous 
calculations, and outward along / = 0. Again, note that the line / = /' = is 
excluded from Region /// by a tubular neighborhood of strictly positive radius, 
so there is an X2 > xi where the integral curve exits Region /// through / = 0. 

Now, consider a point {0,c',X2) along this integral curve with c' > 0. In this 
case, the flow moves such a point rightward. On the other hand, the left boundary 
of Region I moves leftward, approaching / = 0. So there must be an X3 > X2 
such that the integral curve through (0, c', X2) enters the Region I. Collecting our 
findings, we see that every point in C has an integral curve which passes to Region 
/, and therefore corresponds to a solution which is unbounded, and blows up for 
some finite x. □ 

Theorem 55. Suppose (j){x) > and < for all x > 0. The set Z of initial 
conditions to (14.31) that lead to bounded solutions 

1. lies within A = RiU R2 and is 

2. nonempty, 

3. closed, 

4- unbounded, 

5. connected, and 

6. simply connected. 

7. Additionally, the portion of Z corresponding to solutions that enter the inte- 
rior of R2 is a 1-dimensional submanifold of {{f , f , x)\x = 0}. 

Proof. 1. From Lemma [541 all bounded solutions must lie in A. 

2. That there exist bounded solutions in A is the content of Lemmas 152] and 1531 
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3. Now, put Aq = An {{f, f, 0)} and Bq = dA - Aq. Observe that from the 
proofs of the previous theorems, the flow of V along Aq is inward, and the 
flow along Bq is outward. Since the last component of V does not vanish, 
the flow of V causes each point of Bq to lie on an integral curve starting 
on ^0- This establishes a homeomorphism fl from Bq into a subset of Aq. 
In particular, Q is an open map. Now every solution passing through Bq 
is of course unbounded, so Z = Aq — Q{Bq) is evidently closed (it is the 
complement of an open set). 

4. Bq clearly has the topology of Mx [0, oo), so 7ri(-Bo) = 0. Hence, 7ri(f2(i?o)) = 
also, but notice that fl{BQ) contains OAq. Suppose Z were a bounded set. 
Then it is contained in some disk D. But dD is homotopic to a loop in 
Aq — Z, which either lies in int(74o — Z) (in which case the homotopy need 
not move it) or in 9^40. But this means that the loop encloses all of Z, 
and so cannot be contractible in D,{Bq), which contradicts the triviality of 
7ri(n(i?o))- Hence Z is unbounded. 

5. We first show that the portion of Z lying in the region R2 satisfies the 
horizontal line test. First, note that a solution starting in Z n R2 cannot 
exit i?2- For one, it cannot enter Ri, since i?i is an antifunnel. Secondly, 
it cannot exit into — {Ri U R2) since solutions there are all nonglobal. 
Suppose that /i(0) > /2(0) > and /((O) = /^(O) with (/i(0), /((O)) and 
1/2(0), /^(O)) both in Z n i?2. But then 

j^u'M-m) = n{x)-n{x) 

= fl{x)-fl{x)>Q, 

with equality only if /i(0) = /2(0). Hence, ^(/i(a;) - f2{x))>Q for x > 0, 
again with equality only if /i(0) = /2(0). Now all solutions which remain in 
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i?2 are monotonic decreasing and bounded from below, so they must have 
hmits. On the other hand, the only possible limit is (O,lim^_,oo a/ (t>ix)), so 
therefore all bounded solutions in R2 must have a common limit. Therefore, 
we must have that /i(0) = /2(0). Now this means that the portion of Z in the 
region R2 can be realized as the graph of a function from the /' coordinate to 
the / coordinate. Therefore, if Z were not connected, at least one component 
of Z would be a bounded subset, which is a contradiction. 

6. Finally, if Z were not simply connected, the Jordan curve theorem gives 
that there are two (or more) path components to ^^(-Bo) = ^0 ~ Z, which 
contradicts the continuity of Q. 

7. By the connectedness of Z and the horizontal line test in R2, the function 
from the /' coordinate to the / coordinate whose graph is Z (1 int R2 must 
be continuous. Additionally, by the connectedness of Z and the uniquenss of 
solutions to ODE, this implies that the rest of Z whose solutions enter the 
interior of R2 is also a 1-manifold. 

□ 

Definition 56. It is convenient to define, in addition to the initial condition set Z, 
other sets Z^^ C {(/,/', = Xq} such that any integral curve passing through 
a point in Z^jy exists for all x > 0. Similarly, one can define Z'^^. 

Remark 57. If — > as x — 00, we conjecture that Z acquires the structure 
of a 1-manifold with boundary. The series solution (14.111) is not valid at such 
a boundary of Z, since such a solution must remain in Ri and therefore decays 
quicker than the leading coefficient of (14. lip . Indeed, by analogy with the case 
where = 0, the leading term /o of the series solution would vanish, and the 
solution is then asymptotic to — (f){s)ds dt. 
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All solutions in the form of the series solution (I4.1ip enter i?2, so a result of 
this theorem is that one of the two parameters c? or in the series solution is 
superfluous. Since d parametrizes solutions when = 0, we conventionally take 
K = Using this, fl4.17p indicates that a good approximation (as xq ^ oo, locally 
near / = /' = 0) to the set Z^^ is the set 

wfj') = = {u\n\\f = \n- 

Remark 58. If0^P>Oasx— >■ oo, then it is not true that Z is a 1-manifold 
(with boundary). Indeed, Z has the structure of a 1-manifold attached to the 
teardrop-shaped set M from Lemma [Ml 

4.6 Geometric properties of the initial condition set Z 

Lemma 59. Suppose 0(x) > 0, (f)'{x) < for all x > and (p as x ^ oo. 
Then the set Z intersects {(/,/', a;)|/' = 0}. 

Proof. First, observe that Z intersects the boundary of i?i in x = 0, since we have 
by Lemmas [52] and [53] solutions entirely within Ri and its complement. Using 
the fact that Z is connected and the Jordan curve theorem, Z must intersect the 
boundary of Ri in the plane x = 0. This reasoning also applies for each Z^g with 
Xo > 0, so that we can find points in the intersections Z^^ fl dRi for each Xq > 0. 
Also note that for the backwards flow associated to our equation (ie. the flow of 
—V), solutions which enter Ri must exit through the plane x = 0. Hence there 
exists a sequence of points {Fn} C Z with F„ = (/„, /4, 0) such that the integral 
curve through Fn passes through Gn = {gn,gny^) & Zn H dRi for each integer 
n > 0. 
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Discern three cases: 



1. If any Fn are in Quadrants I or II, then since Z is connected, it must intersect 
{f = 0}. 

2. If any F„ are in Quadrant III, observe that the flow across the surface S = 
{(/,/', a;) 1 = |/'^,/' < O} is right-to- left. Thus the integral curve must 
cross into Quadrant II on its way to Gn- Therefore, the set Z cannot intersect 
the surface S, and so it must intersect {/' = 0}. 

3. Assume all the F„ lie in Quadrant IV. Observe that {Fn} is a closed subset 
of -Ri n {x = 0}, which is compact. Hence some subsequence of {Fn} must 
have a limit, say F. Since Z is closed, F E Z. But in the portion of Ri lying 
in the x = plane and in Quadrant IV, we have that 

d 



and 



f=f-^<0 
ax 



ax 



Hence fn > g'n- But since — > 0, (7^ — > 0, so F lies on {/' = 0}. 



□ 

Lemma 60. Under the same hypotheses as Lemma lSR Z also intersects the half 
plane {/ = 0,/'>0}. 

Proof. Using Lemma [59], we form a sequence {Fn} C Z such that the integral 
curve through Fn passes through {/' = 0, f > 0,x = n} for each integer n. (This 
can be done without loss of generality, because if any integral curves pass through 
{/' = 0, / < 0}, then the proof is complete by connectedness of Z.) Note that this 
sequence is entirely contained within Ri by Lemma [SI 
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Discern three cases: 



1. There exists an F„ in either of Quadrants II or III. The result follows by the 
connectedness of Z. 

2. There exists F„ in Quadrant IV. This cannot occur unless the integral curve 
through Fn passes through Quadrant III since the flow along {/' = 0} points 
inward into the portion of Quadrant IV inside R\. 

3. Otherwise, we assume {Fn} is entirely contained within Quadrant I. In this 
case, note that 

-f/' = /^-<^<0. 
ax 

Hence the /'-coordinate of the integral curve through each is positive on 
the interior of Quadrant I. Hence 

;f/ = /'>0, 
ax 

so jn < 9n- But g'n — > since ^ 0, so any limit point of {i^n} will have 
/-coordinate equal to zero. By the compactness of R\V\{x — 0} and the 
closedness of Z, this implies that Z intersects {/ = 0, /' > 0}. 



□ 

Lemma 61. Suppose (f){x) > for all x > 0, (p ^ as x ^ oo, and that there 
exists an Xq > such that for all x > Xq, (f)'{x) < 0. Then the set Z intersects 
{/ = 0,/'>0}. 

Proof. We follow the pattern of proving the existence of an intersection for an 
open interval in x containing xq, and then constructing an a priori estimate for 
the /'-coordinate of this intersection. 
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Apply Lemma [601 to xq, we have that Z^^ intersects {/ = 0, /' > 0}. Let 
(0, /o,a;o) he in this intersection. Note that 

d 



and 



,-/ = /> 
ax 



;f/' = f -0=-0<O 
ax 



when evaluated there. As a result, the integral curve passing through (0, /Q,a;o) 
must pass through Quadrant II first, say for x G (xi, xq). Then evidently, Z^^ must 
intersect {/ = 0,/' > 0}. 

Now since > between xi and Xq, and [a;i,a;o] is compact, there is an 
open set in M.^ containing the intersection of each Z^ with {/ = 0,/'>0} for each 
X G [a;i,a;o], such that in this open set £f <K < 0. As a resuh, /( > Hence 
the /'-coordinate of the intersection point of Z^ with {/ = 0, /' > 0} is decreasing 
with increasing x. (Since we have — </) > — </>, it is decreasing at a rate no faster 
than 0. This implies that this intersection point has /'-coordinate no larger than 
J^" (f){x)dx + /^ at a; = 0.) Now since solutions through Z^^ exist for all a; > by 
definition, this suffices to show that Z intersects {/ = 0, /' > 0}. □ 



Remark 62. The line of reasoning used in the third case of each of Lemmas 
and [60] (and also in [6T1) fails if we try to continue Z much farther. This is due to 
the nonmonotonicity of df'/ dx in Quadrants II and III. More delicate control of (p 
must be exercised to say more. 

Calculation 63. Towards the end of the more delicate results mentioned in Re- 
mark [621 it is useful to know the maximum speed along integral curves on points 
in the region Ri in the /- and /'-directions. By this we mean to compute for fixed 
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X the maximum values of 

I /'I for the /-direction 

(4.18) 

\P — (p{x)\ for the /'-direction 

in The first is easy to maximize: we simply look for the maximum value of /' 
in Ri, which is a maximum of 



/'=\/^/^-2/0(x) + ^03/2(3,) 



for — 2a/0(s) < / < a/ (pix). This occurs at / = — and has the value of 
^8730^/^. For the second part of dHHD, it is easy to see that the maximum is 
30(a;). In summary, 

If I < A/I^^^'^fa;) for the f-direction 



1/ — (f){x)\ < 3(f){x) for the /'-direction 



on Ri. 



Using this calculation, we can impose a stronger bound on the decay of 
and constrain the set Z further. 

Lemma 64. Suppose (j){x) > 0, (f)'{x) < — for all x > for some 
< k < 1 and D > 1. Then the set Z is contained within {/ > — A;a/0(O)} and 
intersects each vertical and horizontal line in {/ > 0} exactly once, and intersects 
{/' = 0} only once. 

Proof. That Z intersects {/ = 0, /' > 0} and {/' = 0, / > 0} at all follows from 
Lemmas [59] and EH Now consider the region A G Ri shown in Figure 14.71 and 
defined by 

A = R,n ({/' > 0, / < u {/' < 0, 2f < 3/'^} 
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Figure 4.7: The region A of Lemma [Ml 

The boundary segments strictly to the right of the boundary labelled 1 in Figure 
14. 71 are evidently inflow, so long as > 0. The boundary labelled as 1 in the figure 
moves with speed 

which is greater than maximum speed in the /-direction given in fl4.19p . This 
implies that the boundary moves faster than any solution inside Ri. Hence it is 
an infiow portion of the boundary. On the other hand, the curved segment of the 
boundary to the left has been shown to be outfiow, in Lemma [521 

We observe that the boundary marked 2 in Figure 14.71 moves with speed 

which is strictly faster than the boundary marked 1 in Figure [4?71 and the boundary 
marked 3 in Figure 14.71 moves with speed 
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noting that f'{-^/^f(x)) = ±^/8/3(l)^/^{x) is the value of the maximum /'- 
coordinate of Ri at a given x value. This last speed is greater than the maximum 
speed in the /'-direction given by fl4.19p since < -Dy/6(l)^/\x). (Notice that 
VQ< since < A; < 1.) 

Since D > 1, this means that both the boundaries marked 2 and 3 in Figure 
14. 71 overtake any solution constrained to be within Ri. As a result, every solution 
within the region A must leave it within finite x. But the only way to leave A 
causes a solution to enter — {Ri U R2), so every solution which contains a point 
in A cannot exist for all a; > by Lemma [511 Therefore, Z is contained within 
U R2) - A. 

Now consider the region B which is defined by 

B = Rin ({/' > 0, / < 0} U {/' < 0, 2f < 3/2}) , 

which is simply the region A, with k taken to be zero. The portion of the boundary 
of B lying in the {/ = 0} plane is inflow. We can therefore apply the reasoning 
of the vertical line test: Suppose (/i, /{, 0), (/2, 0) G Z with /i = /2 > and 
/i > /2- Then we have both (at x = 0) 

d 



and 



^(/l-/2) = /(-/>0. 



which gives that /i — /| > for some open interval about x = 0. Then, -^-(/^ 



f'2) ^ 0, which implies that in fact ^(/i — /2) > 0. However, since all bounded 
solutions tend to the common limit of zero, we have that this implies f[ = / 
at X = 0. (Note that since each solution starts in Z Ci (i?i — B), we have that 
neither solution can become negative, since that would involve entering B G A or 
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leaving RiU R2-) This implies that there is a unique intersection of Z with each 
vertical line. The same reasoning applies in the case of the horizontal line test, as 
in Theorem [551 D 

Lemma 65. Suppose (f){x) > for all x > 0, and that (j)\x) < —D^^(f)^^'^{x) for 
all X > xo > for some < k < 1 and D > 1. Additionally, suppose that for all 
X e [0,Xo], 

xo-x< V^I^^^ , (4.20) 

8 P3/4 
3 



where P = max^gjo^xo] 4>{x)- Then the set Z is contained within {/ > — a/0(O)} and 
intersects each vertical and horizontal line in {/ > 0} exactly once, and intersects 
{/' = 0} only once. 



Proof. The set Z^^ is constrained to lie within the set {/' > — A;a/0(xo)}, by 
Lemma [Ml (replacing xq by zero). Now using the /-direction part of (14.191) . the 
smallest /-value attained in Z^ is 



-(ff^^'^{x)dx — k\/ 0(xo). 



If X < Xo, we have 



<fl\x)dx-k^f^{^ > J p^/\x-xo)-k^/^ 

XQ 1 O V O 



X 



by flCTD . As a result, Z^ C {f > -^^(x)} for each X < xq. This additionally 
means that in the backwards flow, the entire portion of Z^ contained in {/ < 0} 
is moving away from the plane {/' = 0}, which completes the proof. □ 



83 



Remark 66. The condition that (j)'{x) < C(jy'^^{x) imphes 

< C 

~(t)-^'\x) < Cx + C 
(h(x) < 



iC" - x)^ 



for some C" and C" . Notice that this condition is satisfied when the series solution 
converges by Theorem HTl 



4.7 Solutions on the entire real line 

We now combine the results for (14.31) and (14.41) to discuss properties of the solutions 
to (14.21) . When (f){x) is monotonically decreasing for x > 0, we have by Lemma 
[5l] that the initial condition set stays within Ri U i?2- In particular, Z C {/' < 
|0'^''^(O)}. If we relax the restriction of monotonicity, we obtain a similar result. 

Lemma 67. If f = f{x) is a bounded solution to the initial value problem (14.31) 
with G n L°°(R) then /'(O) < ^^8^ 



3/4 
oo • 



Proof. Since / is a solution to (14.31) . then it must satisfy 

/" = f- m > f - 



Now Lemma [Ml shows that all bounded solutions to g" = — ||0||oo lie in the 
closure of the set M given by 

M = !^{9,g')\l9' - \g" - gUU + >o,g< 



Since this set M is bounded, we can find the maximum value of /', which is 

fL. = v^ll0ll^'- □ 
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Lemma 68. Consider solutions to (14. 2p on the real line, with cf) G fl 
If for some — oo < A < B < oo, 

-B 



{x)dx>^^\^i sup + ( sup 

then no bounded solutions exist. 



A V 3 \ x&{-oo,A] x-e[-B,oo) 



Proof. For a solution /, we have that /" = — (pix) > — Integrating both 
sides we have 

f{B)-f{A)>- [%{x)dx. 



A 



By Lemma \67\ bounded solutions on 



• x> B have f'{B) < y |(sup^g(_^^^] and 

• onx> A, they have f'{A) < (sup^gj^,^) |</)(a;)|)^/'^, 

so a necessary condition for there to be a bounded solution is that 
(p{x)dx < \M { sup + ( sup 

A V 3 \ a;e{-oo,A] xe[B,oo) 



□ 



Corollary 69. A necessary condition for bounded solutions to (14.21) to exist if 
(peC^n L°°(R) IS (p{x)dx > 0. 

Proof. Suppose bounded solutions exist. By the proof of Lemma [68| if we let 

/n 
(f){x)dx, 

and 

hn = \[l({ sup + ( sup 



3 



oo,— n] a;£[n,oo) 
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then gn < hn for each positive integer n. But the continuity of hmits gives 



Definition 70. A function G C^nL°°{M.) will be called M-s/ia]5e(i if there exists 
an xq > such that for all |x| > Xq, > and 

• is monotonic increasing for x < —Xq and 

• (f) is monotonic decreasing for x > Xq. 

Theorem 71. Suppose (f) is a positive M-shaped function, then solutions exist to 



Proof. Observe that by Lemma [611 we have that the set Z intersects {/ = 0, /' > 
0}. Additionally, by Theorem EHl we have that Z also lies in R2, which is un- 
bounded in Quadrant IV. Likewise, the set Z' (for (14.41) ) intersects {/ = 0, /' < 0}, 
and becomes unbounded in Quadrant I, so Z fl Z' must be nonempty, and at least 



Theorem 72. Suppose (j) is a positive M-shaped function which additionally sat- 
isfies the decay constraints of Lemma {6^ for x > and x < separately, then a 
unique positive solution exists to (14.21) . (Note that for x < 0, the inequalities and 
signs in Lemma\6R must he reversed, mutatis mutandis.) 

Proof. By the Theorem [7T| there exist solutions to (14.21) . one of which comes from 
the intersection oi Z H Z' in the half-plane {x = 0, / > 0}. The vertical-line test 
in Lemma [65] allows one to conclude that the solution which passes through that 




□ 



one point in this intersection is in the half-plane {x = 0, / > 0}. 



□ 
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-0.25 ' ' ' ' ' ' ' 1 

-0.1 0.1 0.2 0.3 0.4 0.5 0.6 

f 

Figure 4.8: The sets Z and Z' in Example [731 

half-plane must continue directly to the region R2 of Lemma [531 without crossing 
the plane {/ = 0}. Thus this solution is strictly positive. 

On the other hand, Lemma [65] indicates that Z may lie only in Quadrants I, II, 
and IV, while the set Z' must lie in Quadrants I, III, and IV. On the other hand, 
the vertical- and horizontal-line tests ensure a unique intersection of Z and Z' in 
Quadrants I and IV, so the solution is unique. □ 

Example 73. We examine the family (pd^) = ce~^^/^, which is M-shaped when 
c > 0. Notice that when c < 0, then the necessary condition of Corollary [69] is not 
met, so solutions do not exist for all x G M. When c = 0, then the trivial solution 
/ = is the only solution. For c > 0, we examine (p'^x) = —xce"^^^"^. Figure [T8] 
shows the sets Z and Z' for the case when c = 0.05. In particular, one notes that 
there appears to be a unique point of intersection. 
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We find the xq for wliicli larger x satisfy (j)'{x) < —A^/2(|)^^^{x)/{k^/3)^. 



which occurs if z > l^c^/"^, so we may take xq = f^c^/*^. 

By way of example, if we fix xq = 4/3, then k = \/6c^^^. (We enforce < < 1 
by taking c small.) Now we must check to see if fl4.20p holds. In this case, we need 
to see if c can be chosen so that xq — x = 4/3 — xis bounded above by 

p3/4 ~ v^c3/4 



> 



which can be made as large as one likes by taking c sufficiently small. Noting that 
this last line is a constant in x completes the bound. Therefore, there is a unique 
positive solution for = /" — + ce~^'^/^ with c G [0, e) for some e > 0. 

Remark 74. Taken together, the results of Corollary EHl and Theorems [7T] and 
72] for M-shaped provide the following story about solutions to the equation 
= /" — + on the real line: 

• If the portion of (f) where it is allowed to be negative is sufficiently negative, 
then no solutions exist, 

• If is positive, then a solution will exist. There is no particular reason to 
believe that this solution will be strictly positive or unique. 

• If the decay in the monotonic portions of is fast enough, there is exactly 
one solution, which is strictly positive. 
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4.8 Numerical examination 



4.8.1 Computational framework 

Notice that the results of Remark [71] are not sharp: nothing is said if (j) has a 
portion which is negative, but still satisfies the necessary condition of Corollary 
l69l Further, if is positive, but does not satisfy the decay rate conditions, nothing 
is said about the number of global solutions that exist. Answers to these questions 
can be obtained by combining the asymptotic information we have collected about 
the sets Z and Z' with a numerical solver. In particular, we can obtain information 
about the number of global solutions to (14. 2p for any M-shaped 0. 

Implicit in the use of a numerical solver is the following Conjecture: 
Conjecture 75. There are only finitely many smooth global solutions to (14.21) . 

Suppose that is an M-shaped function, and that Xq is such that 0(x) is 
monotonic decreasing for all x > xq and is monotonic increasing for all x < — xq. 
(If decreases fast enough, we can choose Xq so that the series solution converges 
on the complement of {—xq^xq) for sufficiently small initial conditions.) Then we 
have the sets Z'_,j.^^ and Z^^ of initial conditions to ensure existence of solutions 
on (— oo, — Xo] and [xq, oo) respectively. Then any solution to the boundary value 
problem 

f 

= /"(x) — P{x) + 0(x) for — Xo < X < Xq 
< (/(-Xo),f(-Xo))GZi,.„ (4-21) 
(/(xo),/'(xo))GZ.„ 

V 

extends to a global solution of (14.21) . So all one must do is solve (I4.2ip numerically. 
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An easy way to do this is to numerically extend the sets and Z^^ to Z' and 
Z respectively (ie. extend them to the plane x = 0) and compute Z' ^Z. 

In order to analyze (14. 2p numerically, it is necessary to make a choice of 0. 
Evidently, the numerical results for that particular choice of cannot be expected 
to apply in general. However, a good choice of will suggest features in the 
solutions that are common to a larger class of 0. We shall use 

0(x; c) = (x' - c)e-"'/^ (4.22) 

where c is taken to be a fixed parameter. (See Figure 149]) This choice of has the 
following features which make for interesting behavior in solutions to (14. 2p : 

• 0(x; c) > for c < 0. In this case, there are solutions to (14. 2p . by Theorem 
[TTl On the other hand, the decay rate conditions are not met over all of M so 
the uniqueness result of Theorem [72] does not apply. Inded, the decay rate 
conditions are met only for sufficiently large |a;|, but not for |a;| small. 

• If c > is large enough, it should happen that no solutions to (14.21) exist, 
since the necessary condition of Corollary [69] is not met. Indeed, the integral 
of vanishes when c = 1. 

4.8.2 Bifurcations in the global solutions 

Once computed, the numerical solutions can then be tabulated conveniently in a 
bifurcation diagram. That is, consider the set in given by (c, /(O), /'(O)) for 
each solution /. Evidently, by existence and uniqueness for ordinary differential 
equations, each solution can be uniquely represented by such a point. The results 
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Figure 4.9: The function c) for various c values 

of such a computation are shown in Figure I4.10[ In this diagram, the solutions 
are color-coded by the number of positive eigenvalues of ^ — 2/ as an operator 
C^(M) —>■ C"^(M), which will be shown in Chapter [7| to be the dimension of the 
unstable manifold of an equilibrium solution /. (It should be noted that the green 
curve continues for c < —1.2, but was stopped for display reasons.) 

Considering the bifurcation diagram, it appears to indicate that (14. 3 p undergoes 
a saddle node bifurcation at approximately c = 0.7706, and a subcritical pitchfork 
bifurcation at c = 0.0501. The results agree with Theorem [7l|, in that solutions 
do exist when c < 0. The saddle node bifurcation was anticipated by the general 
shape of 0. For c > 0.7706, global solutions do not exist, which was qualitatively 
predicted by Corollary [691 

However, there are some stranger features of the bifurcation diagram. Most 
prominently, the bifurcation diagram appears simply to end near c = —0.4652, 
and at each branch of the pitchfork at c = 0.0740. It is important to verify that 
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Figure 4.10: Bifurcation diagram, coded by spectrum of ^ — 2/: green = 
nonpositive spectrum, blue = one positive eigenvalue, red = two 
positive eigenvalues 

these are not numerical or discretization errors. If these ends are to be thought of as 
valid bifurcations, very likely, ^ — 2/ acquires a zero eigenvalue there. Plotting 
the smallest magnitude eigenvalue gives some credence to this possibility. (See 
Figure 021) 

As another check, one can measure the size of the existence interval for solutions 
to fl4.2p . centered at x = 0. Looking in the (c, /(0))-plane (taking /'(O) = 0), one 
can find the first x such that the solution exceeds a particular value. This is shown 
in Figure I4.13[ in which one sees the same general shape as in the bifurcation 
diagram. (The jagged nature of the graph along the actual bifurcation diagram 
is due to aliasing.) However, for c < —0.4652, the lower branch clearly continues 
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Typical global solution curves 




X 



Figure 4.11: Typical global solutions: green are from the positive branch, 
the blue one is taken from the lower branch with /'(O) = 0, and 
the red ones are from the fork arms 



Smallest magnitude of eigenvalue of numerical approx to (A-2f) 

0.35 I : 1 : 1 : 1 




c 



Figure 4.12: Smallest-magnitude eigenvalue measured along the lower branch 
with /'(O) = 
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Bifurcation diagram overview 




Figure 4.13: Estimate of existence interval length 

into solutions that exist for only finite x. So the end bifurcation indicates a failure 
of the solutions to (14.21) to exist for all x. 

From the point of view of (14.11) (the parabolic problem), the end bifurcations 
indicate that the equilibria are degenerate in the sense of Morse. It is easy to con- 
struct a 1-parameter family of 2-dimensional flows for which end bifurcations occur. 
The resulting equilibrium solutions in that case always acquire a center manifold 
- essentially a zero eigenvalue as noted above. In an infinite-dimensional fiow, 
however, equilibria can be degenerate without having a center manifold. This is a 
manifestation of the fact that infinite-dimensional spaces are not locally compact. 
In Chapter El we show that in fact all equilibrium solutions are asymptotically un- 
stable. However, in Chapter [71 we find that all equilibria have finite-dimensional 
unstable manifolds whose dimension is determined by the dimension of the positive 
eigenspace of ^ — 2/. In particular, there are equilibria whose unstable manifolds 
are empty, yet they are unstable. 
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4.9 Conclusions 



In this chapter, an approach for counting and approximating global solutions to 
a nonlinear, nonautonomous differential equation was described that combines 
asymptotic and numerical information. The asymptotic information alone is 
enough to give necessary and sufficient (but not sharp) conditions for solutions 
to exist, and provides a fairly weak uniqueness condition. More importantly, the 
asymptotic approximation can be used to supply enough information to pose a 
boundary value problem on a bounded interval containing a smaller interval where 
asymptotic approximation is not valid. This boundary value problem is well-suited 
for numerical examination, and the combined approach yields much more detailed 
results than either method alone. 

The techniques and results of this chapter should apply to more general kinds 
of differential equations. Indeed, there should be no particular obstruction to 
extending any of the analysis to equations of the form 



where N > 2. Regions Ri and R2 are then relatively easy to construct, and similar 
results hold for them as are shown here. Additionally, the asymptotic series for 
large x can be obtained since = /" — has easily-found explicit solutions. 

There is considerably more difficulty in trying to understand solutions to equa- 
tions like 



where G is some polynomial with G{0) = 0. In this more general setting, explicit 
solutions to = /" + G{f) are significantly harder to find and work with. There is 
also no reason to expect that Theorem [38] will hold, since G may have several zeros. 



O = r(a:)-r(x)+0(x), 




(4.23) 
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As a result, the asymptotic analysis becomes essentially unavailable. Therefore, 
the only remaining tool is Wazewski's antifunnel theorem, for which one still needs 
a good description of the asymptotic behavior of solutions. 
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Chapter 5 

Existence of nontrivial eternal solutions 
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5.1 Introduction 

The existence of eternal solutions poses a potentially difficult problem, because the 
backward-time Cauchy problem is well known to be ill-posed. Obviously, equilib- 
rium solutions are trivial examples of such eternal solutions, and in Chapter H] we 
showed that they exist. It is not at all clear that there are other eternal solutions, 
and indeed there may not be. In this chapter we use a pair of nonintersecting equi- 
hbrium solutions to construct a heteroclinic orbit which connects them. Therefore, 
the set of heteroclinic orbits is generally nonempty. 

As has been done in previous chapters, we will work with the more limited 
equation (11. 6p 



where is a certain smooth function which decays to zero. In particular, when 
there are at least two equilibria whose difference is never zero, there exist nonequi- 
librium eternal solutions to (15.11) . By an eternal solution, we mean a classical 
solution that is defined for all t G M and x G M. We follow the general tech- 
nique for constructing "ancient solutions," which was used in a different context 
by Perelman. 

5.2 Equilibrium solutions 

We choose (f){x) = (x^ — 0.4)e-^ /2. It has been shown in Chapter H] (see Figure 
I4.10p . that in this situation, there exists a pair of equilibrium solutions /+, /_ with 
the following properties: 
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1. /+ and /_ are smooth and bounded, 

2. and /_ have bounded first and second derivatives, 

3. /+ and /_ are asymptotic to Q/x^ for large x, and so both belong to L^(]R), 

4. > /-(x) for all 

5. there is no equilibrium solution /2 with > f2{x) > f-{x) for all x, 
and additionally, there exists a one-parameter family Qc of solutions to 



2. go = /_ and 0o = 0, 

3. 4>a{x) < 4>b{x) and ga{x) > gb{x) for all x ii a > h. 

The latter set of properties can occur as a consequence of the specific structure 
of /_. For instance, consider the following result. 

Proposition 76. Suppose /_ G C^'"(M) satisfies the above conditions and addi- 
tionally, there is a compact K (ZM. with nonempty interior such that /_ is negative 
on the interior of K and is nonnegative on the complement of K . Then such a 
family gc above exists. 



Proof. (Sketch) Work in T/-_C^'"(]R), the tangent space at /_. Then (15 ■2p becomes 
its linearization (for he, say), namely 



= g'cix) -9l{x) +(t)c{x) 



(5.2) 



with 



1. ce [0,1), 



= /l','(x)-2/_(x)/le(x) + (</., -(/.). 



(5.3) 
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Consider the slighly different problem, 



= y"{x) - 2f_{x)y{x) + v{x)y{x), 



(5.4) 



where f is a smooth function to be determined. If we can find a f < such that 
y > and y ^ as \x\ ^ oo, then we are done, because we simply let vy = (f>c — (p 
in (15. 3p . In that case, he = y has the required properties. We sketch why such a v 
exists: 

• If V = 0, then y = is a solution, giving = /_ as a base case. 

• If v{x) = —2\\u\\ool3{x) for (3 is a. smooth bump function with compact sup- 
port and l3\K = 1, then the Sturm-Liouville comparison theorem implies 
that y has no sign changes. We can take y strictly positive. However, in this 
case, the Sturm-Liouville theorem imples that ther are no critical points of 
y either, so y may not tend to zero as |x| oo. 

• Hence there should exist an s with < s < 2||-u||oo such that if v{x) = 
—sl3{x), then y has no sign changes, one critical point, and tends to zero 
as |x| —>■ oo. This choice of v is what is required. (The precise details of 
this argument fall under standard Sturm-Liouville theory, which are omitted 
here.) 



In what follows, we shall not be concerned with the exact form of 0, but rather 
we shall assume that the above properties of the equilibria hold. Many other 
choices of (p will allow a similar construction. 

Lemma 77. The set 



□ 



W = {ve C^{R)\f-{x) < v{x) < U{x) for all x} 



(5.5) 
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is a forward invariant set for ( 15. ip . That is, if u is a solution to (IS.ip andu{tQ) G 
W, then u{t) G W for all t > to. 



Proof. We show that the flow of fl5.ll) is inward whenever a timeslice is tangent to 
either /_ or /+. To this end, define the set B 

B = {v E C^(]R)|/_(x) < v{x) < for all x, and there exists an xq 

such that f (xq) = f+{x) or f (xq) = 

Without loss of generality, consider a. v E B with a single point of tangency, 
'^{.^o) = f-i^o)- At such a point Xq, the smoothness of v and /_ implies that 
Av{xo) > A/_(xo) using the maximum principle. Then, if m is a solution to (15.11) 
with u{0,x) = v{x), we have that 

^^^^ = Av{xo)~v'{xo)+<f>{xo) 

> A/_(xo)-/!(xo) + 0(xo) = O, 



hence the flow is inward. One can repeat the above argument for each point of 
tangency, and for tangency with /+ as well. □ 

Lemma 78. Solutions to the Cauchy problem 

'-^ = '^-u\t,x)+m, 
M(0,a;) = U{x) e Wc 

where 



(5.6) 



W^ = {ve C\R)\gc{x) < v{x) < /+(x) for all x} 

for c G [0, 1) have the property that they lie in L} fl L°°(R) for all t > 0. We shall 
assume that U has hounded first and second derivatives. 
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Additionally, when c G (0,1), solutions to fl5.6p cannot have /_ as a limit as 
t — > oo. 

Proof. The fact that solutions he in nL°°(M) is immediate from Lemma 1771 and 
the asymptotic behavior of /+, /_ (Section l4.4p . Observe that for each c G [0, 1), 
Wc is forward invariant, and that Wa C Wb if a > 6. Since /_ is not in Wc for c 
strictly larger than 0, the proof is completed. □ 

The following is an outline for the rest of the chapter. All solutions to (15.61) 
have bounded first and second spatial derivatives. This implies that all of their first 
partial derivatives are bounded (the time derivative is controlled by (15. ip ). Using 
the fact that (15.11) is autonomous in time, time translations of solutions are also 
solutions. We therefore construct a sequence of solutions {uk} to Cauchy problems 
started at t = 0, Ti, T2, ... which tend to /+ as t ^ +00, but their initial conditions 
tend to /_ as A; ^ 00. By Ascoli's theorem, this sequence converges uniformly on 
compact subsets to a continuous eternal solution. 

5.3 Integral equation formulation 

In order to estimate the derivatives of a solution to (15. 6p . it is more convenient to 
work with an integral equation formulation of (15.61) . This is obtained in the usual 
way. 
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Am - + 
-v^ + (t) 

(|-a)"w-.^) 

/oo nt poo 

H{t,x -y)U{y)dy+ / H{t - s,x - y) {(t>{y) - u'^{s,y)) dyds, 
oo Jo J —oo 

(5.7) 

where H(t, x) = ^=e 4t is the usual heat kernel. 

Calculation 79. We begin by estimating the first derivative of u for a short time. 
Let T > be given, and consider < t < T. The key fact is that / H(t, x)dx = 1 
for all t. Using ([5TD 



d_ 

di 



du 
'di 

A 1 u 



u 



du 
dx 



< 



< 



< 



< 



dU 




dx 


oo -^0 


dU 




dx 


oo JQ 


dU 




dx 


oo ^0 


dU 

dx 


+ T 

oo 



f f°° d 



'0 J-oo 

t POO 



d_ 

dy 

Hit - s,x-y) 



dx 



{H{t-s,x-y)){(P{y)-u\s,y)) 

dy ds 

uy J 

2||'u||f^ / -— ds. 



dy ds 



d(j) 2 ^ 
dy dy) 





du 


f 

Jo 


dx 



This integral equation fence is easily solved to give 



du 




dU 


+ T 







^ ( 




dx 


oo \ 


dx 


oo 


dx 



,2tmax{|l/+|U,|l/_|U} 



K2T 



Calculation 80. With the same choice of T as above, we find a bound for the 
second derivative in the same way: 
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< 



< 





+ T 

oo 




dx"^ 


dx"^ 


d'U 


+ T 




dx^ 


dx"^ 



+ 

oo ^0 



d / du 
dy V dy 



ds 



du 



dx 



+ 2\\u\ 



dx"^ 



ds 



for some which depends on U, 0, and T. 

Calculation 81. Now, we extend Calculation [80] to handle t > T, 



dx'^ 



< 



dx^ 



+ 
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00 

f't poo 

/ / H{t-s,x-y){(f){y)-u'^{s,y))dyds 

Jt J-oo 

H(t — s, x] 
ds + K'^ 



dx^ 



T poo 



H{t - s,x- y){(t){y) - u^{s,y))dyds 



J -00 



+ 



dx'^ 



dx"^ 
* 1 



S\ s 



m\i + \\A\i)ds 

ds 



X S \/ s 



hence there is a uniform upper bound on 



which depends only on the initial 



conditions, 0, and T. 



Lemma 82. Let f G C2(R) be a bounded function with a bounded second deriva- 
tive. Then the first derivative of f is also bounded, and the bound depends only on 
. and II /"II 00 ■ 



Proof. The proof is elementary. The key fact is that at its maxima and minima, 
/ has a horizontal tangent. From a horizontal tangent, the quickest /' can grow 
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is at a rate of ||/"||oo- However, since / is bounded, there is a maximum amount 
that this growth of /' can accrue. Indeed, a sharp estimate is 



ll/lloo < V2II/IUII/II00. 

□ 

Using the fact that u is bounded. Lemma [82] imphes that the first spatial 
derivative of u is bounded. By (15. ip . it is clear that the first time derivative of u 
is also bounded. 



Lemma 83. ^45 an immediate consequence of Lemmas 78 and{8^ the action in- 
tegral 

1 du 1 
J -00 2 ^ 

is bounded. Therefore, the solutions to the Cauchy problem (15. 6p all tend to limits 
as t —s> 00 ( Corollary By Lemma \78l we conclude that they all tend to the 

common limit of /+ when c > 0. 

Proof. The latter two terms are bounded due to the fact that u lies in Li fl 
for all t. The bound on the first term comes from combining the fact that u and 
its first two spatial derivatives are bounded with the asymptotic decay of f±, and 
is otherwise straightforward (use L'Hopital's rule). □ 



5.4 Construction of an eternal solution 



Let 

Uk{x) = (1 - 2-''~')gyk{x) + 2-'-'f+{x), for > 



105 



noting that Uk f- as k ^ oo. Since Uk is a convex combination of /+ and 
gi/k, it follows that Uk G Wi/k for all k. Also, since /+ and /_ have bounded first 
and second derivatives, the {Uk} have a common bound for their first and second 
derivatives. 

Now consider solutions to the following set of Cauchy problems 



Uk{Tk,x) = Uk{x). 



(5. 



We choose Tk so that for all A; > 0, ma;(0,0) = uo(0,0). We can do this using 
the continuity of the solution and Lemma [83l As — >■ cxd, solutions are started 
nearer and nearer to the equilibrium /_, so we are forced to choose — oo as 

fc — >• CXD. 

It's clear that each solution Uk is defined for only t > Tk- However, for each 
compact set S C M^, there are infinitely many elements of {uk} which are defined 
on it. The results of the previous section imply that {uk} is a bounded, equicon- 
tinous family. As a result, Ascoli's theorem implies that {uk} converges uniformly 
on compact subsets to a continous u, which is an eternal solution to (15.11) . 

Our constructed eternal solution will have the value m(0, 0) = mo(0, 0), which is 
strictly between /+ and /_. As a result, the eternal solution we have constructed 
is not an equilibrium solution. By Lemma [831 it is a finite energy solution, so it 
must be a heteroclinic orbit connecting /_ to /+. 
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Chapter 6 

Instability of equilibria 
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6.1 Introduction 



(This chapter is available on the arXiv as [36].) 

If we try to apply standard Morse theory to our semilinear parabolic equation, 
we encounter a serious difficulty. In particular, the stability of the linearization 
about an equilibrium of our system is not sufficient to ensure that the equilibrium 
is stable, even though there may not be a zero eigenvalue. The instability of 
equilibria for (11.21) is not a new fact, having been studied carefully in the 1980s. 
This chapter is included for completeness, providing an explicit construction of 
a sequence of solutions starting near the equilibrium which all blow up. Indeed, 
there is a complementary result of stability in certain weighted norms, described 
in [15] and [38]. 

Note that the right side of (II. 2p is an operator which has a spectrum which 
includes zero, so stability is possible (as in the unforced heat equation), though 
not guaranteed. This is in stark contrast to the situation in finite-dimensional 
settings, where asymptotic stability of the linearized system implies stability of 
the equilibrium (in particular, zero is not an eigenvalue of the linearized operator). 
(See [5], for instance.) Essentially, this difficulty suggests that each critical point 
is degenerate. (This line of reasoning is completed in Chapter [71) This implies 
that Morse theory (even when strengthened to its natural infinite-dimensional form 
[31] ) cannot be used to study the dynamics of our system. 

Again, we study the general situation by working with the simpler Cauchy 
problem (II. 6p . (The computations we exhibit in this chapter will carry over mutatis 
mutandis to (II. 2p .) However, it is useful to center on an equilibrium solution. That 
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is, we apply the change of variables u{t, x) ^ u{t, x) — f{x) to obtain 
§iu{t, x) = ^u{t, x) - 2f{x)u{t, x) - u'^{t, x) 
< u{0,x) = h{x) G C°°(R) (6-1) 
t > 0,x G M, 

where / G C^(R) is a positive function with two bounded derivatives. (By C^, 
we mean the space of smooth functions which decay to zero.) We interpret / as 
being an equilibrium of the original problem (11. 6p . 

The assumption that / be positive requires some motivation. With this assump- 
tion, the zero function is an asymptotically stable equilibrium for the linearized 
problem, 

^u{t, x) = -^uit, x) - 2f{x)u{t, x). (6.2) 

by a standard comparison principle argument. This corresponds neatly to the case 
of zero positive eigenvalues of ^ — 2/ which was found numerically in Figure 
14.101 Intuition would suggest that this implies / is a stable equilibrium of (II. 6p . 
However, using a technique pioneered by Fujita in [18], we will show that this 
equilibrium is not stable in the nonlinear problem, even if the initial condition has 
small p-norm for every 1 < p < oo. 

Fujita showed that if / = 0, then the zero function is an unstable equilibrium 
of (16. ip . The cause of the instability in (16. ip is the decay of /, for if / = const > 0, 
then the comparison principle shows that the zero function is stable. We extend 
Fujita's result, so that roughly speaking, since / — away from the origin, the 
system is less stable to perturbations away from the origin. Another indication 
that there may be instability lurking (though not conclusive proof) is that the 
decay of / means that the spectrum of the linearized operator on the right side of 
(16. 2p includes zero. 
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6.2 Motivation 



The problem (16. ip describes a reaction-diffusion equation [T2], or a diffusive lo- 
gistic population model with a spatially-varying carrying capacity. The choice of 
/ positive means that the equilibrium u = describes a population saturated at 
its carrying capacity. Without the diffusion term, this situation is well known to 
be stable. The decay condition on / means that the carrying capacity diminishes 
away from the origin. 

The spatial inhomogeneity of / makes the analysis of (16. ip much more compli- 
cated than that of typical reaction-diffusion equations. The existence of additional 
equilibria for (16. ip is a fairly difficult problem, which depends delicately on /. (See 
[6] for a proof of existence of equilibria in a related setting.) 

6.3 Instability of the equilibrium 

Given an e > 0, we will construct an initial condition h G C°°(M) for the problem 
(16. ip . with \\h\\p < e for each 1 < p < oo, such that ||M(t)||oo — * oo as t — 
T < oo. In particular, this implies that m = is not a stable equilibrium of 
(16. ip . at least insofar as classical solutions are concerned. We employ a technique 
of Fujita, which provides sufficient conditions for equations like (16. ip to blow up 
|18j . (Additionally, [TT] contains a more elementary discussion of the technique 
with a similar construction.) Our choice for h can be thought of as a sequence 
of progressively shifted gaussians, and we will demonstrate that though each has 
smaller p-norm than the previous, the solution started at h still blows up. 
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6.3.1 The technique of Fujita 



The technique of Fujita examines the blow-up behavior of nonhnear parabohc 
equations by treating them as ordinary differential equations on a Hilbert space. 
Suppose u(t) solves 

^ = Lu{t) + N{u{t),t), (6.3) 

where L is a linear operator not involving t, and N may be nonlinear and may 
depend on t. Suppose that v(t) solves 

dv{t) 



dt 



-L*v{t), (6.4) 



where L* is the adjoint of L. Let J(t) = {v(t),u(t)). We observe that if \ J(t)\ oo 
then either ||i^(t)|| or ||M(t)|| also does. So if v{t) does not blow up, then we can 
show that ||M(t)|| blows up, and perhaps more is true. If we differentiate J{t), we 
obtain the identity 

pit) = j^{v{t)Mt)) 

= {-L*v{t), u{t)) + {v{t), Lu{t) + N{u{t),t)) 
= {v{t),N{uit),t)), 

where there is typically a technical justification required for the second equality. 
It is often possible to find a bound for {v{t), N{u{t),t)) in terms of J{t). So then 
the method provides a fence (in the sense of for J{t), which we can solve to 
give a bound on |J(t)|. As a result, the blow-up behavior of u(t) is controlled by 
the solution of an ordinary differential equation (for J{t)) and a /mear parabolic 
equation (for v{t)), both of which are much easier to examine than the original 
nonlinear parabolic equation. 
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6.3.2 Instability in W for 1 < j9 < oo 



We begin our application of the method of Fujita by working with L = ^ — 2/ 
and N{u) = —v? in (16.31) . Since (16. 4p is then not well-posed for all t, we must be a 
little more careful than the method initially suggests. For this reason, we consider 
a family of solutions f ^ to (16.41) that have slightly extended domains of definition. 
It will also be important, for technical reasons, to enforce the assumption that the 
first and second derivatives of / are bounded. 

Definition 84. Suppose w = w{t,x) solves 



w(0, x) = Wo{x) > 0. 

Define v^{s, x) = w(t — s + e,x) for fixed t > and s < t + e. Notice that by the 
comparison principle, v^{s,x) > 0. 

Lemma 85. Suppose that w solves (16. 5p . Then w, |^ G Co(M). 

Proof. The standard existence and regularity theorems for linear parabolic equa- 
tions (see m], for example) give that w,^,^ e L^{R) and that w G C^{R). The 
comparison principle, apphed to and §1^^ gives that the first and second 
derivatives of w are bounded for each fixed t. (This uses our assumption that / 
has two bounded derivatives.) 

The lemma follows from a more general result: if (7 G nL^(M) for 1 < p < 00 
and g' G L°°(M), then g G Co(]R). To show this, we suppose the contrary, that 
lima;^oo g{x) 7^ (and possibly doesn't exist). By definition, this implies that there 
is an e > such that for all a; > 0, there is a ?/ satisfying y > x and \g{y)\ > e. Let 
S = {y\ \g{y)\ > e}, which is a union of open intervals, is of finite measure, and 




(6.5) 
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has supS* = oo. Let T = {y\ \g{y)\ > e/2}. Note that T contains S, but since g' is 
bounded, for each x ^ S, there is a neighborhood of x contained in T of measure 
at least e/||(7'||oo- Hence, since supT = sup S = oo, T cannot be of finite measure, 
which contradicts the fact that g G L*'(M) with 1 < p < oo. □ 

Lemma 86. Suppose -u : [0, T) x M — > R zs a classical solution to (16. ip with u < 
and u{t) G L°°(R) for each t G [0, T). Then 

— w(t,x)h{x)dx < ( / -|| — , ... dsj , (6.6) 

where w is defined as in Definition 

Proof. Define 

J.(.) = (6,7) 

First of all, we observe that since u G L°^(R), v^{s, ■)u{s, ■) is in L^(R) for each 
s <t. 

Now suppose we have a sequence of compactly supported smooth func- 

tions with the following properties: [29] 

• m„ G C°°(M), 

• mn{x) > for all x, 

• supp(m„) is contained in the interval (— n — l,n + 1), and 

• m„(x) = 1 for < n. 

Then it follows that 

J^{s) = lim / v^{s, x)u{s, x)mn{x)dx. 
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Now 



— lim / v^{s,x)u{s,x)mn{x)dx 

as n— »oo J 

lim lim — / {v^{s + h,x)u{s + h,x) — Ve{s,x)u{s,x))mn{x)dx. 

h—>0 n— >oo h J 



We'd like to exchange limits using uniform convergence. To do this we show that 
lim lim — / (^£(5 + h, x)u{s + h,x) — Ve{s, x)u{s, x)) mn{x)dx (6.8) 

n— ioo /i— »0 h J 

exists and the inner limit is uniform. We show both together by a little computa- 
tion, using uniform convergence and LDCT: 



lim lim — / (i'e(s + /i, x)u{s + h,x) — Ve{s, x)u{s, x)) mn{x)dx 

n— >oo /i— >0 h J 

[ f d d \ 

— lim / [ —vJs,x)u(s,x) + vJs,x)—u(s,x) ] 'mn(x)dx 
J yds OS J 



lim 



dx"^ 



v^{s, x) + 2/(a;)fg(s, x) ) x)mn{x) + 



d 



v^{s, x) I -Q-^u{s, x) — u (s, x) — 2f{x)u{x) J mn{x)dx 



— lim / —Ve{s, x)u'^{s, x)mn{x)dx. 

n— >oo J 



Minkowski's inequality has that 

VeUrUndx 



< J Ve\u\mndx < v^mndx^ VeU^nindx^ , 



since Ve,mn > 0. This gives that 
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(J v^urrindxY 
J Vefriidx ' 



hence the inner Umit of (16.81) is uniform. On the other hand, 

\v,{s,x)u^{s,x)mnix)\ < v,{s,x)\\u{s)\\l^ E L\R) 

so the double hmit of (16.81) exists by dominated convergence. Thus we have the 
fence 

ds - \\v,is)\\,- ^ ■ ^ 

We solve the fence (16.91) to obtain (note < 0) 

1 ^ dJAs) 1 



||u£(s)||i ds {Je{s)) 

* 1 . 1 1 



2 



-ds < 



\\v,{s)\\^ - Ut) J,(0) 
ds < 



/o \\v.is)h - MO) 
Taking the limit as e ^ of both sides of the inequality yields 



t 1 / rt 



w(t, x)h(x)dx < ( / - — ; TT^'^s I = ( [ T, — , , 11 ds 

^ ^ ^ ^ -\Jo \\wit-s)h J yjo ik(s)iii 

as desired. □ 

Remark 87. Since we are interested in proving the instability of the zero function 
in (16.11) . consider u{0,x) = h{x) = — e for e > 0. Then (16.61) takes on the simple 
form 



e 



'PPJrdsKl. (6.10) 
\Ms)\\i 
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So in particular, \\u{t) ||oo blows up if there exists a T > such that e 




1. 



The stability of the zero function in (16 .ip depends on the stability of the zero 
function in (16.51) - the linearized problem. If the zero function in the linearized 
problem is very strongly attractive, say ~ e^*, then 



and so a small choice of e < 1 does not cause blow-up via a violation of (I6.10p . On 
the other hand, blow-up occurs if it is less attractive, say ~ for a > 0. 

Because then 



whence blow-up occurs before t = 

In the particular case of f{x) = for all x, we note that w is simply a solution 
to the heat equation, which has ||w(t)||i = \\wo\\i for all t (by direct computation 
using the fundamental solution, say), so blow up occurs. Thus we can recover a 
special case of the original blow-up result of Fujita in [18]. 

Theorem 88. Suppose a sufficiently small e > is given. Then for a certain 
choice of initial condition h{x) with \\h\\p < e for all 1 < p < oo, there exists a 
T > for which limt^y- ll^(^)l|oo = oo. 

Proof. First, it suffices to choose ||u(0)||i < e and ||n(0)||oo < e, since 



We assume, contrary to what is to be proven, that ||M(t)||oo does not blow up for 
any finite t. In other words, assume that u : [0, oo) x M ^ M is a classical solution 
to (16. ip . with ||M(t)||oo < oo for all t. We make several definitions: 






i/p 
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Choose < P < mill |e, 

Choose 7 > small enough so that 



K, (6.11) 



2772 

for some some arbitrary K > 1. 

Since < / G C^(M), we can choose an xi such that 

f{x) < 7 when x < Xi. (6-12) 

Next, we choose Xq < xi so that 

V^II/lloo(l-erf(^^)) <7 (6.13) 

for all < t < Notice that any choice less than xq will also work. 
Choose the initial condition for (16. ip to be 

m(0, x) = h{x) = -/?e^'^'(--o)^ (6.14) 
This choice of initial condition has ||u(0)||oo = P < 11^(0) 111 = 2n^/'^P^/^ < 



e, and 



jj, = 2/5^/2. (The value of fi will be important shortly.) 



Finally, let Wo{y) = 6{y — Xq) (the Dirac 5-distribution), and suppose that 
w solves (16. 5p . In other words, choose w to be the fundamental solution to 
(16. 5p concentrated at xq. Note that the maximum principle ensures both 
that w{t,x) > for all t > and x G M and that ||w(t)||i < ||w(0)||i = 1 for 
alH > 0. This allows us to rewrite (16. 6p as 

-t j w{t,x)h{x)dx <1. (6.15) 
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Now we estimate the integral in (]6.15p . Notice that 



— / w{t, x) {—h{x)) dx 



dx'^ 
dx 



■ 2f{x)w{t, x) j {—h{x))dx 
2 + 2f{x)h{x) I w{t^x)dx, 



where Lemma [HS] ehminates the boundary terms. Now suppose z solves the heat 
equation with the same initial condition as w, namely 



dz 

at 



(6.16) 



2;(0, x) = Wo{x) = 6{x — Xq). 



The comparison principle estabilishes that z{t, x) > w{t, x) for alH > and a; G M, 
since f,w>0. As a result, we have that 



d 
dt 



w{t, x) {—h{x)) dx > 



dx^ 



+ 2f{x)h{x) z{t,x)dx 



> -/i-2/3 / f{x)z{t,x)dx, 



where /x 
result, 



1^(0) and P = ||'u(0)||oo, which is an integrable equation. As a 



w{t,x) {—h{x)) dx > (3 — nt — 2(3 I I I f{x)—^=e ^^wo{y)dydxds. 




(6.17) 
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On the other hand using our choice for wq, 

f{x) e 4» Wo{y)dy dx ds = / / f{x)—^^=e 4s rfxc/s 

< y -^=== y e 4» cia;+ ll/llooy e cia; 1 cis 



4 2"-'"°°J„ V, 2v^ 



4 2"-'""i„ V 2v^ 

< ^<7Vi. 



we have used fl6.12p . (16. 131) . and assumed that < t < Then (16.151) becomes 

2^3/2^3/2 



+ Pt, 



V27K 

using our choices of /i, 7, and initial condition. Maple reports that the maximum 
of A{t) = -2[3^/'^t^ - '^'^^l^f^'^ + is unique, occurs at < to < 4^, and has the 



asymptotic expansion 



Thus for all small enough e > /?, we obtain a contradiction to fl6.15p since K > 1. 
Thus, for some T < to < 00, lim^^r- ||'u(t)||oo = 00. □ 



6.4 Discussion 



Theorem [88] gives a fairly strong instability result. No matter how small an initial 
condition to (16. ip is chosen, even with all p-norms chosen small, solutions can blow 
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up so quickly that they fail to exist for all t. This precludes any kind of stability for 
classical solutions. Like the analogous result in Fujita's paper, the kind of initial 
conditions which can be responsible for blow up are of the nicest kind imaginable 
- gaussians in either case! 

It must be understood that the argument in Theorem [88] depends crucially on 
the decay of /. Without it, the lower bound on / w(t,x){—h{x))dx decreases too 
quickly. Indeed, if / = const > and h{x) > — /, then the comparison principle 
demonstrates that the zero function is asymptotically stable. On the other hand, 
any rate of decay for / satisfies the hypotheses of Theorem [88], and so will cause 
()6.ip to exhibit instability. 

Finally, although we have examined the case where the nonlinearity in (16.11) 
is due to M^, there is no obstruction to extending the analysis to any nonlinearity 
like \u\^, with degree k greater than 2. A higher-degree nonlinearity would result 
in a somewhat different form for (16.61) . but this presents no further difficulties to 
the argument. Indeed, by analogy with Fujita's work, higher-degree nonlinearities 
would result in significantly faster blow-up. 
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Chapter 7 

Cell complex structure for the space of 
heteroclines 
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7.1 Introduction 



In this chapter, we determine that all unstable manifolds of (11.61) 

du{t,x) d'^u{t,x) 



dt dx"^ 



(7.1) 



are finite dimensional. This is not a particularly new result, indeed Theorem 
5.2.1 in [20j can easily be made to apply with the Banach spaces we shall choose. 
Theorem 5.2.1 in [20] shows the existence of a smooth finite dimensional unstable 
manifold locally at an equilibrium. One can then use the iterated time-1 map of the 
flow for (17. ip to extend this local manifold to a maximal unstable manifold. There 
are also finite Hausdorff dimensional attractors for the forward Cauchy problem 
on bounded domains [32]. However, we shall exhibit a more global approach to 
the finite dimensionality of the unstable manifolds. This approach allows us to 
examine the finite dimensionality of the space of heteroclinic orbits connecting a 
pair of equilibria, which is a new result in the spirit of [15j. The techniques used 
here depend rather delicately on both the degree of the nonlinearity (quadratic) and 
the spatial dimension (1). Both of these are important in the standard methodology 
as well, as the portion of the spectrum of the linearization in the right half-plane 
needs to be bounded away from zero. In the case of (17.11) . the spectrum in the 
right-half plane is discrete and consists of a finite number of points. 



7.2 The linearization and its kernel 

We begin by considering an equilibrium solution / to (17. ip . As discussed in Chapter 
m this solution has asymptotic behavior which places it in C^flL^ nL°°(M). We are 
particularly interested in solutions which lie in the a-limit set of /, those solutions 
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which are defined for alH < and tend to /. As in previous chapters, center on 
this equihbrium by applying the change of variables u{t,x) u{t,x) — f{x) to 
obtain 

a;) = ^u{t, x) - 2f{x)u{t, x) - u^{t, x) 
u{0,x) = h{x) G C\R) 

(7.2) 

\imt^_^u(t,x) = f{x) 
t < 0,x e M. 

Thus we have a final value problem for our nonlinear equation. All solutions to 
(17. 2p will tend to zero as t ^ — oo uniformly by Theorem [28l Of course, (17.21) is 
ill-posed. We show that there is only a finite dimensional manifold of choices of h 
for which a solution exists. 



7.2.1 Backward time decay 



The decay of solutions to zero is a crucial part of the analysis, as it provides the 
ability to perform Laplace transforms. In the forward time direction, one obtains 
upper bounds for solutions by way of maximum principles, and lower bounds for the 
upper bounds by way of Harnack estimates. In the backward time direction, these 
tools reverse roles. Harnack estimates provide upper bounds, while the maximum 
principle provides lower bounds for the upper bound. In the proof of Theorem [281 
the latter was used to some advantage. In this section, we briefiy apply a standard 
Harnack estimate to obtain an exponentially decaying upper bound. 

Harnack estimates for a very general class of parabolic equations are discussed 
in [28] and [IJ. In those articles, the authors examine positive solutions to 

du 

div A(a;, t, m, Vu) — = B{x, t, u, Vm), 
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where x G M", and A : M^^^^ ^^^^ q . ^2n+2 _^ ^ satisfy 

|A(x,t,M,p)| < a|p| + c|m| + e 

\B{x,t,u,p)\ < h\p\+d\u\ + f 

p-A(x,t,u,p) > -\p\'^ — dlul"^ — g, 
a 



for some a > and b, ...g are measurable functions. For a solution m defined on 
a rectangle i?, the authors define a pair of congruent, disjoint closed rectangles 
R^, C R with i?~ being a backward time translation of R^. The main result 
is the Harnack inequality 



maxM < 7 ( minw + L ] , (7.3) 
R- \R+ J 

where 7 > depends only on geometry and a (but not b, ...g) and L is a linear 
combination of e, /, g whose coefficients depend on geometry. 

In the case of ( 17. 2p . or indeed of the analogous equation with higher degree 
terms, we have that (17. 3p will apply with L = 0. Notice that the conditions on 
A, B are satisfied because any solution to (17. 2p is automatically a finite energy 
solution, and therefore is bounded and has bounded first derivatives. The only 
difficulty is that (17.31) applies for positive solutions, while (17. 2p may have solutions 
with negative portions. However, one can pose the problem for the (weak) solution 
of 



d\u\ 



sgn {u) (Au — — 2/n) 



dt 

= A|m| — u\u\ — 2f\u\ 



> A|m| - -2|/||m 
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for which we only get positive solutions. By iterating (I7.3P we have that solutions 
to (17.21) decay exponentially as t — — oo. 



7.2.2 Topological considerations 

Definition 89. Let Ya{X) be the subspace of ^^(X, CO'°(M)) which consists of 
functions which decay exponentially to zero like e"*, where < a < 1. We define 
the weighted norm 



M|y. = ||e""*lk(t)||c;o,. 



and the space 

Ya{X) = {u = U{t,x) G C^(X,C°'°(M))|||M||y„ < Oo} . 

In a similar way, we can define the weighted Banach space Za{X) as a subspace of 
C°(X,CO'°(R)). It is quite important that Ya and Za are Banach algebras under 
pointwise multiplication. 

In light of the previous section, solutions to (17. 2p are zeros of the densely defined 
nonlinear operator N : Ya{{—oo,0]) oo,0]) given by 

About the zero function, the linearization of N is the densely defined linear map 
L : y„((-oo, 0]) ^ Zaii-oo, 0]) given by 

d d 

where we define H = ^ — 2f. Also note that L is the Frechet derivative of N, 
which follows from the fact that Ya and Za are Banach algebras. 
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Remark 90. We are using C°'"(]R) instead of C°(M) to ensure that and L be 
densely defined. We could use space of continous functions which decay to zero, 
or the space of uniformly continous functions equally well. 

Convention 91. We shall conventionally take a > to be smaller than the small- 
est eigenvalue of H. 

We show two things: that the kernel of L is finite dimensional, and that L 
is surjective. These two facts enable us to use the implicit function theorem to 
conclude that the space of solutions comprising the a-limit set of an equilibrium 
is a finite dimensional submanifold of Y'a((— oo, 0]). 



7.2.3 Dimension of the kernel 

Lemma 92. // / is an equilibrium solution, then the operator L : ya((— oo,0]) 
Za{{— 00,0]) in (17. 5p has a finite dimensional kernel. 

Proof. Notice that the operator L is separable, so we try the usual separation 
h{t,x) = T(t)X{x). Substituting into (17. 5p gives 

f d \ 



T X 



X 



for some A G C. The separated equation for T yields T = C^e'^^. Since we are 
looking for the kernel of L mYa C. L°°(R^), we must conclude that A must have 
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nonnegative real part. On the other hand, the spectrum oi H = — 2fj 
is strictly real, so A > 0. Indeed, there are finitely many positive possibilities 
for A each with finite-dimensional eigenspace. This is a standard fact about the 
Schrodinger operator H since / is an equilibrium (Proposition lllll) . Thus L has 
a finite dimensional kernel. □ 



7.2.4 Surjectivity of the linearization 

In order to show the surjectivity of L, we will construct a map F : Za{{—oo, 0]) — > 
Ya{{— 00,0]) for which L oT = idz^- That is, we construct a right-inverse to L, 
noting of course that L is typically not injective. We shall derive a formula for F 
using the Laplace transform t> h-i> u 

v{s,x) = / e^^v(t,x)dt, 

where 3?(s) > —a and v G oo,0]). 

Since Lemma [92] essentially solves (17. 2p . we will be solving the inhomogeneous 
problem with zero final condition 

^ - ^ + 2f{x)vit, x) = -wit, x) e Z,((-oo, 0]) 

(7.6) 

v{0,x) = 

for t < 0. The Laplace transform of this problem is 

d'^'vis, x) 

svis,x) ^ —4. 2fix)v{s,x) = w(s,x) 

ox^ 

{H + s)v{s,x) = W{s,x). 

Choose a vertical contour C with > 3?(s) > —a, so that the Laplace transforms 
are well-defined, and that the contour remains entirely in the resolvent set of —H. 
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c 



-a 



Spectrum of -H 



Figure 7.1: Definition of the contour C 



Tlien we can invert to obtain 



f (s, x) = (if + s) ^w{s,x). 



Using tlie inversion formula for tlie Laplace transform yields 

v{t,x) = — I e'^'^H + sy^w{s,x)ds 
21X1 Jc 

= — [ e-'\H + s)-^ I e'^w{T,x)dTds 
2vr« Jc Jt 



We can obtain operator convergence of the operator-valued integral in paren- 
theses if we deflect the contour C. Choose instead the portion C of the hyperbola 
(See Figure El]) 

(where A is the smallest magnitude eigenvalue of —H) which lies in the left half- 
plane as our new contour. Then, since —H : C°'" C°'" is sectorial about 
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(A — a)/2 (Proposition 11121) . Theorem 1.3.4 in [20] implies that the integral 

defines an operator- valued semigroup e"^^'^"*), so the formula for T is given by 

T{w){t,x) = e'^^-^-'^w{T,x)dT. (7.8) 

It remains to show that the image of F is in fact Ya, as it is easy to see that its 
image is in L°°. That the image is as advertised is not immediately obvious because 
the contour deflection C ^ C changes the domain of the Laplace transform. In 
particular, the derivation given above is no longer valid with the new contour. 

Therefore, we must estimate ||f (recall that A is the smallest magnitude 
eigenvalue of —H) 



< 



< 



< 



1 

27n 



c 



CO 



1 





K(s+a)t / 5R(s+a)r| 



27r Jc' |s — A| 
Ki\\w\\z, 

27r 
Ki\\w\\z^ 



w\\z dr ds 



■t \s — A| 



-5R(s+a)t_ 



1 _ ^ns^a)i^ 



[s + a 



TT 



ds 



c, \s-\\\^{s + a)\ 



< Ko\\w\ 



where < Ki,K2 < oo are independent of t and w. We have made use of the 
estimate in Proposition 11121 of the norm of {H + s)^^ : C^'" when s is in 

the resolvent set of —H. In particular, note that the choice of C being to the left 
of —a is crucial to the convergence of the integrals. Thus the image of F lies in 
Za- The backward-time decay of ^ is immediate from the Harnack inequality, so 
in fact the image of F lies in Ya. 
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Theorem 93. The linear map L : Ya{{—oo,0]) —>■ oo,0]) is surjective and 

has a finite dimensional kernel. Therefore the set A^^^(O) is a finite dimensional 
manifold, which is the unstable manifold of the equilibrium f . The dimension of 
A^~^(0) is precisely the dimension of the positive eigenspace of H . 

Proof. The only thing which remains to be shown is that the domain Ya sphts into 
a pair of closed complementary subspaces: the kernel of L and its complement. 
That its complement is closed follows immediately from a standard application of 
the Hahn-Banach theorem. (Extend id^er l to all oiYa-) □ 

Combining the fact that an equilibrium solution can have an empty unsta- 
ble manifold (we numerically computed the dimension of the eigenspaces of L in 
Chapter H]) and is yet unstable, we have proven the following result. 

Theorem 94. All equilbrium solutions to (17.11) are degenerate critical points in 
the sense of Morse. 

7.3 Linearization about heteroclinic orbits 

We can extend the technique of the previous section to the linearization about a 
heteroclinic orbit. The resulting generalization of Theorem |93] is that the connect- 
ing manifolds of (17. ip are all finite dimensional. 

Suppose that m is a heteroclinic orbit of (17.11) . Let /-,/+ be the equilibrium 
solutions of (11.21) to which u converges as t ^ — oo and t +oo respectively. 

Suppose that Aq : M ^ (0, oo) is the smallest positive eigenvalue of H{t). It 
is easy to see that Aq is piecewise C^, for instance, see Proposition 1.7.2 in [25]. 
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Propostion IllOl ensures that Aq is a bounded function. We will define a pair of 
bounded, piecewise functions Ai and A2 which will aid us in defining a two more 
pairs of function spaces. Let Ai : R ^ (0, 00) be a bounded, piecewise function 
with bounded derivative which has the following properties: 

• Ai(t) is never an eigenvalue of H{t), 

• limt^_oo ^ < 1, and 

• since u —>■ f± uniformly, for a sufficiently large i? > 0, Ai can be chosen so 
that there are no jumps on its restriction to M — [—R, R]. 

Defining A2 is a somewhat more delicate problem. We would like to exclude 
the solutions which lie in the unstable manifold of /+, since they cannot lie in 
the space of heteroclines from /_ —>■ f+- We do this by separating the eigenvalues 
corresponding to the intersection of the unstable manifolds of /_ and /+ from those 
which lie in the stable manifold of /+. However, there is an obstruction to this 
technique. In particular, the eigenvalues of H{t) = — 2u{t) vary with time, 
and can bifurcate. To avoid this issue, we need some kind of regularity for the 
eigenvalues to prevent them from bifurcating. We follow Floer [14] in the following 
way: 

Conjecture 95. There is a generic subset (a Baire subset) of choices for the 
coefficients in (11.21) so that if m is a heteroclinic orbit, all of the eigenvalues of 
H{t) are simple. 

Numerical evidence, as exhibited in Chapters H] and [8] suggests that the above 
Conjecture is true. When we assume that all of the eigenvalues of H{t) are simple. 
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- t 



Figure 7.2: Definition of Ai and A2 



and tfierefore do not undergo any bifurcations otlier tlian passing tlirougli zero, we 
sliall say u is a lieterocline contained in Ureg- 

Let A2 be in ^^(R) sucli tliat 

• A2 = Ai on [R, 00), and 

• A2(t) is not an eigenvalue of H{t) for any t. 



We can do this when u G Ureg- See Figure [721 

Definition 96. Define the Banach algebra YxX^) (for i = 1,2) to be the set of u 

in C^(X,CO'"(R)) such that the norm 



\U 



< 00, 



where X is an interval containing zero. Likewise, we can define the spaces Zx^ {X) C 
C*'(X, C°'"(M)) in a similar way. That these are Banach spaces follows from the 
boundedness of the Aj. It is also elementary to see that these are Banach algebras. 
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We then consider Ni,Li as VAi(I^) —>■ Zx^{W), where Li is the hnearization of 
Ni about u for i = 1,2. (Again, since Yx- and Zx^ are Banach algebras, Lj is the 
Frechet derivative of Ni.) For a i G {1,2}, consider the restriction L~ of Li to a 
map YAi((— 00,0]) Zxi{{—oo,0]). We rewrite 

d 92 



Likewise, we can define L^ : yAi([0, 00)) ^Ai([0, C)o)). 

We define the positive eigenspaces for the equihbria as well 

V^if±) = span jt; G C°'"(M)| there is a A > with - 2/±^ = At;| . 

(7.10) 

Note in particular that dim V^{f±) < 00. 

Lemma 97. If u E U^eg is a heterocline that converges to f± as t ^ ±cxd, then 
the operator Li has a finite dimensional kernel for i G {1,2}, and in particular 

lim dim V'^iuit)) — hm dim V^{u{t)) < dim ker Li < dim ker L^ < 00. 

t—*~oo t— ►+00 

(The condition u G Ureg is only necessary for the i = 2 case.) 

Proof. Notice that the first term of (17.91) has finite dimensional kernel by Lemma 
[92] and closed image by Theorem [931 The second term of (17.9!) is a compact 
operator since m ^ /_ uniformly. Thus L~ has a finite dimensional kernel. Let 
span{t>m}m=i = kerL~ and consider the set of Cauchy problems 

§ = 1^ - 2m/i for t > 

at 0.^ ^^^^^^ 
h{0,x) = Vm{0,x). 

Standard parabolic theory gives uniqueness of solutions to (17.111) . and that a so- 
lution h lies in the kernel of Lf, the restriction of Li to [0, 00) x M. Therefore 
dim ker L,- < dim ker L~ < 00. 
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For the other inequahty, modify u outside of [—R, R]xM. to get a m so that the 
hnearization Li of N about u satisfies 

• ker Li is isomorphic to ker Li as vector spaces, 

• m|(-oo -i?)xR = /-, and 

• lf|(iJ,oo)xK = /+• 

We can do this for a sufficiently large R, since u tends uniformly to equilibria. 
Then the flow of 

dh d'^h 
9^=9^ + '^^ 

defines an injective linear map from the timeslice at —R to the timeslice at R. 
(That is, it gives an injective map from C'^'°'(M) to itself - injectivity being an 
expression of the uniqueness of solutions.) Each element v of the kernel of Lj 
evidently must have f (— -R) G and v{R) ^ \^+(/+). Therefore, the injec- 

tivity ensures that the intersection of the image under the fiow of V^{f^) with the 
complement of has at least dimension dim V^+(/_) — dim V^{f\). □ 

Remark 98. Multiplication by u, C^(R2, C°'"(]R)) ^ ^"(M^) is not a compact 
operator, in particular note that dim ker L'l = oo. 

Theorem 99. Let u be a heterocline of (17.11) which connects equilibria f±. There 
exists a union |JM„ of finite dimensional submanifolds o/ C"'^(M, C°'°(]R)) 
which 

• contains u and 

• consists of heteroclines connecting /_ to /+. 

G Ureg, then Mu has dimension lim(_^„oo dim ^"''(^(t)) — lim(_^oo dim V^{u{t)), 
and this is maximal among such submanifolds M^. 
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Proof. Observe that Li is surjective, since it is easy to show that the formula 



converges, where we note that the contour changes with time. As it happens, the 
computation in [2U] goes through with the only change that at t = 0, we deflect 
the contour to the right, rather than the left (as in Figure [7!T]) . Since Lemma [971 
shows that Li has finite dimensional kernel, then it follows that = N^^{0) is 
a union of finite dimensional manifolds, with a finite maximal dimension. It is 
obvious that consists entirely of heteroclinic orbits and contains u. 

It remains to show that the dimension of is as advertised and maximal. 
Observe that L2 is a compact perturbation of an operator L'2 : 5^A2(I^) Zx^iM) 
which is time-translation invariant. This follows from the precise choice of A2 being 
continous and not intersecting the eigenvalues of H. L2 and are both surjective 
by exactly the same reasoning as for Li. is injective by using separation of 
variables as in Lemma [22] (noting that all nontrivial solutions blow up in the 1^2 
norm). Therefore the Fredholm index of L^, hence L2 is zero. However, this implies 
that L2 is injective. 

Since L2 is bijective, any solution to L2U = which decays faster than e/^2(t)(it 
as t ^ —00 ends up growing faster than '^'^^^'''^^ as t — +00, and in particular 
does not tend to zero. As a result, such a solution cannot be in ker Li. This 
implies that dim ker Li < limf^_oo dim V^{u{t)) — limf^oo dim V~^{u{t)), which 
with the estimate in Lemma [HZI completes the proof. □ 

Remark 100. Even if m ^ Ureg (when there exist nonsimple eigenvalues of H{t)), 




is a well defined right inverse of Li. This involves showing that 
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the function Ai can still be constructed. As a result, we always get that the 
connecting manifold is finite-dimensional. 

Corollary 101. The space of heteroclinic orbits has the structure of a cell com- 
plex with finite dimensional cells. This cell complex structure is evidently finite 
dimensional if there exist only finitely many equilibria for ( I7.ip . 



7.4 Conclusions 



We have shown that the tangent space at an equilibrium splits into a finite dimen- 
sional unstable subspace, and infinite dimensional center and stable subspaces. 
However, it is quite clear by Chapter [6] that the center subspace is nonempty and 
large. Indeed, considering the work of [38j, the center and stable subspaces are not 
closed complements of each other. Additionally, we have given conditions for the 
space of heteroclinic orbits to have a finite dimensional cell complex structure. 
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Chapter 8 

An extended example 
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8.1 Introduction 



Consider the following equation 

where the choice of in (11 .Op has been fixed. The bifurcation diagram for the 
equilibria of (18.11) can be found in Figure 14.101 Based on the Theorem [931 the 
number of positive eigenvalues shown in Figure 14.101 corresponds exactly to the 
dimension of the unstable manifold of each equilibrium. 



8.2 Frontier of the stable manifold 

According to Figure H.lOl when c = —1.2, there is only one equilibrium, /q. It 
has empty unstable manifold, though of course it is asymptotically unstable (as 
is shown in Chapter [6]). On the other hand, /o has an infinite dimensional stable 
manifold, which is not all of C°'"(M), as a consequence of the asymptotic instability. 
As a result, its stable manifold has a frontier in C°'"(]R) (which may not be a 
boundary in the sense of a manifold with boundary). We are interested in the 
qualitative behavior of solutions near and along this frontier. We know by Theorem 
[28] that if they tend to /o uniformly on compact subsets, then they do so uniformly. 
It is enlightening to use a numerical procedure to this end. We start solutions at 
the following family of initial conditions 

UA{x) = fo{x) + Ae-^'/''. (8.2) 

Using the Fujita technique (exactly as shown in Chapter [6]), we can show that for 
sufficiently negative A, the solution started at ua will not be eternal. As a result. 
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Figure 8.1: Behavior of solutions near the frontier of the stable manifold of 
/o (horizontal axis is x) 

the family of initial conditions ua intersects the frontier of the stable manifold of 
/q. An approximation to the value of A which corresponds to the frontier can be 
easily found using a binary search. Some typical such solutions are shown in Figure 
18. H and the approximate value of A corresponding to the frontier is A ?a —2.15 

The qualitative behavior shown in Figure 18.11 indicates that there is some kind 
of traveling disturbance in the frontier solutions, which seems like a traveling wave. 
However, such a solution also appears to tend uniformly on compact subsets to /o, 
so in fact it converges uniformly. (The uniform convergence is not obvious from 
the figure, due to the numerical solution being truncated at a finite time.) The 
leading edge of this disturbance collapses to — oo in finite time for solutions just 
outside the stable manifold of /q. 

8.3 Flow near equilibria with two-dimensional unstable 
manifolds 

Also of interest is the structure of the flow in the unstable manifold of the "fork 
arms" which occur at c = 0.0740, as they approach the pitchfork bifurcation at 
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Figure 8.2: Flow in the unstable manifold of a "fork arm." c = 0.0600 (left); 
c = 0.0501 (right) 



Evolution of a heteroclinic orbit Evolution of tfie spectrum of H(t) 




t 



Figure 8.3: A typical heteroclinic orbit to the left of boundary A, with the 
spectrum of H{t) as a function of t. 



c = 0.0501. Figure 18^ shows a schematic of the flow based on numerical evidence. 
Of particular interest is the behavior near the boundary marked A. Solutions to 
the right of the boundary are not eternal solutions - they fail to exist for all t. 
Solutions to the left of A are heteroclinic orbits connecting the equilibrium with 
an unstable manifold of dimension 2 to the equilibrium with an unstable manifold 
of dimension zero. A typical such solution is shown in Figure 18.31 
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Eigenfunction Eigenfunction 




Figure 8.4: Eigenf unctions describing unstable directions at /i 

To examine solutions near the boundary A, we center our attention on the 
case c = 0, which has two equilibria, one of which (call it /i) has a 2- dimensional 
unstable manifold. (This corresponds to the right pane of Figure 18. 2[ ) If we 
hnearize about /i, the operator H = ^ - 2fi : CO'°(M) CO'"(R) has a pair of 
simple eigenvalues, as is easily seen in the right pane of Figure [831 at t = 0. One of 
these eigenvalues is smaller, to which is associated the eigenfunction ei in Figure 
18.41 The eigenfunction 62 is associated to the larger eigenvalue. In Figure 18. 2^ ei 
corresponds to the horizontal direction, and 62 corresponds to the vertical direction. 
From the proof of Lemma 1921 it is clear that {61,62} spans the tangent space of 
the unstable manifold at /i. Therefore, we specify initial conditions UA,e{x) for a 
numerical solver using 

UA,e{x) = fi{x) + A {ei{x) cos 6' + 62(0;) sin 6') . (8.3) 
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(Taking A small allows us to approximate solutions which tend to /i in back- 
wards time.) Since the perturbations along 61,62 are quite small, and indeed the 
eigenvalue associated to 61 is much smaller than that associated to 62, examining 
the numerical results of evolving UA,e is quite difficult. The behavior along the 
boundary occurs at a much smaller scale than /i, yet is crucial in determining 
the long-time behavior of the solution. To remedy this, the boundary behavior 
is better emphasized by plotting UA.e{t,x) — fi{x) instead. Figure [831 shows the 
results of evolving initial conditions (18.31) for A = 0.1 and various values of 6. 

Solutions in Figure 18.51 show a similar kind of behavior as in the case of the 
frontier of /q. There is a traveling front, which moves very slowly in the negative 
x-direction. However, the behavior is quite a bit more delicate. The determining 
factor in locating the frontier of /o is the perturbation in a direction roughly like 
62, which has a large eigenvalue. On the other hand, for /i. Figure 18^ indicates 
that such a direction is not parallel to the boundary of the connecting manifold. 
(The boundary direction is some linear combination of 61 and 62, with a numerical 
value for the angle 6 being roughly 1.114975 radians.) The eigenvalue associated 
to 61 is roughly ten times smaller, and therefore perturbations in that direction are 
much more sensitive. Additionally, the action of the flow is therefore primarily in 
the direction of 61, which tends to mask effects in other directions. For this reason, 
it was visually necessary to postprocess the numerical solutions by subtracting /i 
from them. Otherwise the presence of the traveling front was unclear. 
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9=1.11494 0=1.11496 9=1.11497 




t t t 



9=1.11498 9=1.11499 9=1.115 




t t t 



Figure 8.5: Difference between equilibrium fi and the numerical solution 
started at ua,$, where black indicates a value of -0.2, and white 
indicates 0.2. The horizontal axis represents t, and the vertical 
axis represents x. A = 0.1 in all figures. Starting from the upper 
left, e = 1.11494, 1.11496, 1.11497, 1.11498, 1.11499, 1.115. 
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Chapter 9 

Conjectures and future work 
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9.1 Conjectures about the present problem 
9.1.1 Analytical conjectures 

It seems that under reasonable conditions on the coefficients of fll.2p . all eternal 
solutions ought to be heteroclinic orbits. An easy calculation with the formula (11.31) 
for the action A{u{t)) shows that if |A(-u(t))| blows up, then one of the following 
is true: 

1. \\uit)\\i ^ oo, 

2. ||M(t)||oo — > OO, or 

3. \\Du{u) II2 ^ oo. 

Essentially, eternal solutions which are not heteroclinic orbits are big in some 
sense. Of course, traveling fronts satisfy the ffist condition. On the other hand, 
the Harnack inequality seems to imply that eternal solutions do not blow up in 
the 00-norm as t — * —00. More intriguingly, [45j and [38j show that under certain 
conditions on the coefficients of (11. 2p . global solutions to the forward Cauchy 
problem have a universal bound on their 00-norm. However, these results are 
obtained under the hypothesis that the solution u is strictly negative, a condition 
that is essential to their analysis. Relaxing this condition leads to currently open 
problems. 

Conjecture 102. Suppose all of the coefficients in (II. 2p decay sufficiently fast 
as |x| 00. Then all eternal solutions are bounded in the 00-norm by a universal 
bound, which depends only on the Oj. 
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More ambitious is the following (which involves proving a universal bound for 
the 1-norm as well): 

Conjecture 103. Suppose all of the coefficients in (11 ■2p decay sufficiently fast 
as oo. Then all eternal solutions are heteroclinic orbits. 

Related to both of these conjectures is the conjecture that under suitable decay 
conditions on the a^, there exist only finitely many equilibria (Conjecture [75]) . 

9.1.2 Conjectures related to the topology of the space of 
heteroclinic orbits 

Much of what remains to be understood about the space of heteroclinic orbits of 
(ll.2p and (11.61) involves a more precise understanding of the gluing maps between 
the cells in its cell complex structure. The eventual goal is to construct a homology 
theory, called a Floer homology, for the space of heteroclinic orbits. This would 
allow the space of heteroclinic orbits to be decomposed as a complex of connecting 
manifolds (without boundary) and boundary maps which associate higher dimen- 
sional manifolds to lower dimensional ones. From the outset, degeneracy in the 
sense of Morse provides the biggest obstacle to this kind of theory. In particular, 
nondegeneracy allows one to show that generically, connecting manifolds can only 
have boundaries of one dimension lower. However, in the example of the previous 
chapter, namely that of (II. 6p with = (x^ — cje'^'^^'^, such a statement is still 
true. Perhaps it is possible that one can find conditions for connecting manifolds 
to have codimension-1 boundaries, even in the face of degeneracy in the equilibria. 
Or put another way. 
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Conjecture 104. When the flow of ( ll.2p is restricted to the space of heterochnic 
orbits Ti, all of the equilibria become nondegenerate critical points in the sense of 
Morse. 

Another obstacle is that there needs to be some kind of compactness result for 
the space of heteroclinic orbits modulo time translation (or perhaps modulo action 
of the flow). In Floer's case, he was able to employ Gromov's compactness results 
for pseudoholomorphic curves. This leads to his "no bubbling theorem" . However, 
no such result is known in our case. The closest available results are those of [15] 
and [38], which only hold for positive solutions to (11. 2p that have been centered 
on an equilibrium. It is easy to show that if there exists a positive equilibrium for 
(II. 6p . then their results suffice to show compactness, but the situation of general 
sign is currently an open problem. 

To summarize, we have the following conjectures: 

Conjecture 105. The space of heteroclinic orbits of (11.21) modulo time translation 
is compact in Yx. 

Conjecture 106. There is a generic subset (a Baire subset) of choices for the 
coefficients in (11.21) so that if m is a heteroclinic orbit, all of the eigenvalues of 
H{t) are simple. 

We then define Ck{R) to be the free i?-module generated by the /c-dimensional 
connecting manifolds of (II. 2p . Probably it is best to think of = Z/2, at least to 
fix ideas. 

Conjecture 107. For a generic subset of coefficients in (11.21) . there is a collection 
of maps dk : Ck ^ C'fc-i such that 



147 



• dk is an i?- module homomorphism, for each > 0, 

• 5o = 0, 

• dk-i o dk = for each A; > 0, and 

• V = dk{u) if and only if roughly speaking f is a sum of boundary elements 
of M, obtained by a deformation retraction of a neighborhood of f in m onto 
V. In Floer's case, as is likely in ours, this is a rather involved construction 
called the "gluing theorem" described in [T6] . 

This would turn {C^:{R), 9*) into what is most reasonably called a "Floer com- 
plex." One can then define the Floer homology modules, Fk = ker 5fc/im dk+i and 
formulate the following (reasonable) conjecture: 

Conjecture 108. Let E be the space of heteroclinic orbits of fll.2l) . For a generic 
subset of ttj in (11.21) . Fk = Hk{E;R). That is, the Floer complex computes the 
homology of the space of heteroclinic orbits. 

9.2 Future work on related problems 

9.2.1 Higher spatial dimensions, with decay conditions en- 
forced 

Of course, the most obvious dependence on 1-dimensional space is the equilib- 
rium analysis of Chapter |H The analogous nonlinear elliptic problem (11.71) is not 
well understood. Indeed, very little is known about (I1.7P at all, especially if the 
solutions are allowed to be of general sign. 
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One thing is likely: the spatial decay of heteroclinic orbits is much slower - 
not in L^(R"). Worse, Sturm-Liouville theory is no longer available to control the 
eigenvalues oi H = A — 2u. Therefore, there might be infinitely many positive 
eigenvalues of the operator if, which accumulate at zero. As a result, equilibria 
might have infinite dimensional unstable manifolds. The analysis of the structure 
of the connecting manifolds will therefore not work, though there should be a 
filtration structure based on Lyapunov exponent for the unstable manifolds, which 
should allow for an infinite dimensional cell complex with finite dimensional cells. 
Additionally, the slower spatial decay will disrupt the finite energy classification 
scheme in Chapter [31 

9.2.2 Relaxation of decay conditions on the coefficients 

If we no longer require that the coefficients decay to zero as |x| oo, then (11.21) 
can support traveling wave solutions. Indeed, there can be extremely complicated 
and delicate traveling wave structures if the spatial dimension is also greater than 
1. This will remove the uniform convergence to equilibria, of course. Also, likely 
is that what will be found is that the space of heteroclinic orbits is an infinite- 
dimensional cell complex, perhaps were the "cells" are Banach manifolds. The 
resulting dynamics can therefore be expected to become extremely complicated. 
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Appendix A 

Spectrum of Schrodinger operators 
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A.l Introduction 



This appendix recounts a few standard facts about the structure of the spectrum 
of the Laplacian and Schrodinger operators A, i/ : C^(]R) C°(]R) respectively. 
Nothing in this appendix is original, but it is useful to have the facts and the 
requisite calculations available for reference. 

A. 1.1 Spectrum of the Laplacian operator 

Proposition 109. The spectrum of A : C^(R) C°(R) is the closure of the 

negative real axis. 

Proof. The spectrum of A contains all A e C for which (A — A) is not injective. 
In other words, it contains the solutions to the equation 

u" -Xu^O. 

An elementary calculation yields that u — c\e^^ + C2e~^^. If A is real and non- 
positive, there are nontrivial bounded solutions (which are oscillatory or constant). 
Otherwise the nontrivial solutions are unbounded. Hence the spectrum must 
contain the closed negative real axis. 

Next, we show (A — A)"-*^ exists and is bounded away from the closed negative 
real axis. To show that (A — A)" exists, we find an inversion formula, which is 
valid when A does not he on the closed negative real axis. To this end one can 
solve 

u — Xu = f 
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using a slightly modified version of Calculation [8], 



u{x) 



(-l)e-^(f-^) if ^{\^){y - x) > 
g-v^(y-a;) otherwise 



dy. 



This inversion formula defines a bounded inverse for (A — A) when A is not on 
the closed negative real axis. To see this, simply observe that 



is independent of x. This fact implies that < -^aII/IIc" for some finite K by 

differentiation under the integral. Therefore the complement of the closed negative 



A. 2 Spectrum of Schrodinger operators 

The previous section can be generalized to the case of Schrodinger operators H = 
[A — V) : C^(]R) —>■ C°(M) to obtain a few results of interest. Assume that \^ is a 
smooth function which satisfies lim|2;|^oo V{x) = A. 

Proposition 110. The spectrum of H contains the portion of the real axis less 
than or equal to —A. All of the eigenvalues of H are contained in the portion of 
the real axis less than or equal to the supremum of V. 

Proof. Of course, the eigenvalues A are those where there are nontrivial solutions 
to the equation 




real axis is in the complement of the spectrum of A. 



□ 



u"{x) - {V{x) + X)u{x) = 0. 



(A.l) 
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Recast f[0) first-order system, namely 

A / \ 



d 

dx 



( \ ( 

u 



u 



Tlx) 



( \ 

u 



(A.2) 



1 

Observe that the eigenvalues of T{x) are purely imaginary if A is real and less 
than — V[x). Thus, the fiow restricted to a plane of constant x consists of periodic 
orbits if and only if A < — V[x) (ignoring the origin, of course). 

Now suppose that A is real and A < —A. Then there exists an i? > such that 
for all |x| > R, \ f{x) - A\<\\\^ A\. So on M - \-R, i?], solutions to (TO) will 
all tend to limiting cycles. On [— i?, -R], solutions grow exponentially fast, at a rate 



of no more than a/II^ ||oo + |A|. Thus there exist nontrivial bounded solutions to 
( 1A.1I1 . which are obviously in . Hence the spectrum of H contains (— cxd, 



On the other hand, if A is real and A < sup V < X, the origin is always a saddle 
point for the first order system flA.2p . Thus all solutions to flA.ip are unbounded. 
Likewise, if A G C — M, then the system flA.2l) is a spiral, so [H — A) is injective 
for A ^ (— oo,sup V]. □ 



It is best to treat the portion of the spectrum lying between —A and sup V 
using Sturm-Liouville theory. Indeed, ( lA.lj) on [0, oo) with boundary conditions 

n(0) = a, lim u{x) = 

x— >oo 

is a classic Sturm-Liouville problem. It is known that the eigenvalues of this 
problem are discrete and accumulate only at zero. 

Proposition 111. Suppose V G C°(]R) is positive outside a compact interval. 
Then the operator H has finitely many eigenvalues greater than —A. 
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Proof. It suffices to show that any solution u to 

u -{V + A)u = 

has finitely many zeros. Notice that (assuming by hypothesis that A > —A) 

-V-X<-V-A 

so the Sturm-Liouville comparision theorem states that any eigenfunction v of H 
with eigenvalue A has strictly fewer zeros than solutions to the A = —A case. Also, 
the zeros of solutions which are bounded for half intervals (—00, xq) are monotonic 
in A. 

By hypothesis, V is positive on M — [— -R, R] for some R > 0. By comparision 
with the case of V = on M — [—R,R], zeros of u only occur on [—R,R]. By 
comparison with ||V^||oo, the number of zeros is proportional to -Ra/|| V^||oo, and 
therefore finite. □ 

It is Proposition II 111 that ensures that the finite dimensionality results of Chap- 
ter [7| hold. In the case of higher spatial dimensions, Sturm-Liouville theory does 
not apply (at least if there is no assumed symmetries in the equilibria). It is there- 
fore possible that there is no finite dimensionality for the space of heteroclines for 
higher spatial dimensions. 

It will be technically important in Chapter [7| that {H + s) = {A — {V + s)) is 
sectorial about any real s not in the spectrum of H. It is then useful to consider 
the densely defined operator {H + s) : CO'"(M) ^ CO."(M) intead of ^ C^, 
where < a < 1. 

Proposition 112. If s G M not in the spectrum of H = {A — V), then {H + s) : 
CO."(M) CO'"(M) zs sectorial forO <a<l. 
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Proof. We have already shown that the structure of the spectrum is favorable for 
the sectoriality of {H + s) = (A — (V^ + s)), and the operator is densely defined if 
a > 0. (It is not densely defined if a = 0.) What remains is that the operator norm 
of the inverse must decay like K/\s\, and that its image must consist of continous 
functions. For the former: 



u" — {V + s)u 
u 



uix] 



f 

Vu + f 



2v^ 



(Vu + f) 
-\'-y\^\V{y)u{y) + f{y))dy, 



where s is not in the spectrum. Using Calculation [9l 



\u\ 



< — 

|s| 

J_ 

1^1 



-1 



A - / (Vu + f) 



< 



1- nll^l 



Since V and / are assumed to decay to zero and the kernel e~'^~^'^ is in L^, it 
is immediate that u must also decay to zero. The following calculation shows that 
image of {A — (V + s))~^ consists of Lipschitz functions. Assume Xi < Xq, so that 
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\u{xo) - u{xi] 



2^s 



/XQ poo 
e(y-^o)V-s(yu + f)dy- / e^'^'-y^^'iVu + f)dy 
-oo J xo 

/XI poo 
^{y-xi)V-s^Vu + f)dy+ / e(^i-^)^(yii + /)dy 
■oo J Xl 



2v^ 



+ 
+ 



XI 

oo 



{Vu + f)dy 



Jxo 

+ / e^'''-y^V-s(Vu + f)dy 

J XI 



< 



2v^ 



|^^ + /||c 



p{y-xo)Vs _ ^{y-xi)^/s 



dy 



+ 



Xo 



^{xi-y)Vs _ ^{xo-y)Vs 



dy ^ ^il - e^"°^^') 



< 



Vu + fl 



2^s 



r,Xl^ 



+ 



^_(l_g(xi-xo)7i-) 



< —\\Vu + f\\ 
\s\ 



g(xi-xo)v/i _ 

as |xo — xi| ^ 0. 



But since we've chosen Xi < xq, the above calculation proves that u is Lipscliitz. 
Looking at the Lipschitz constant and the bound on u, it is immediate that the 
C°'^-operator norm of {H + s)~^ decays hke Thus {H + s) is sectorial □ 



It is important to remark that the above proof shows that in fact {A — {V + s))^^ 
is a bounded operator C°(M) C^(M). However, the norm of the operator does 
not decay as ^{s) — > oo as required for a sectorial operator. In particular, 



du 
dx 



< 



i^l||^^ + /||c 
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and 



<K2\\Vu+fl 



for Ki,K2 independent of s and /. Examples can be constructed to show that 
these bounds are tight. 
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